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fully  developed,  turbulent  channel  flow.  Results  show  that  the  mean-shear  (MS) 
and  turbulence-turbulence  (TT)  partial  pressures  and  tt^^,  respectively)  are 
the  same  order  of  magnitude.  The  buffer  region  dominates  most  of  the  wavenumber 
range;  the  viscous  shear  layer  is  significant  at  the  highest-wavenumbers;  buffer 
and  logarithmic  regions  are  important  at  low-wavenumbers.  The  dominant  length 
scales  of  the  MS  source  term  indicate  that  it  is  generated  by  near-wall  shear  layers. 
Over  most  of  the  wavenumber  range,  the  contribution  from  the  buffer  region  is  the 
dominant  TT  component;  in  the  low-wavenumber  range  the  viscous  shear-layer, 
buffer  region,  and  logarithmic  region  are  dominant;  in  the  medium-  and  high- 
wavenumbers  the  viscous  shear-layer  and  buffer  region.  The  most  important  TT 
partial  pressures  are  and  from  the  buffer  region.  It  is  hypothesized 

that  and  are  generated  by  quasi-streamwise  vortices  that  are  parallel  and 
tilted  with  respect  to  the  wall,  respectively,  rnay  be  due  to  near- wall  shear 
layers  and  spanwise  vortices,  but  is  much  less  important  than 
and  Tx'^^  from  the  viscous  shear-layer  are  the  dominant  high  wavenumber  partial 
pressures;  they  may  be  due  to  the  downward  side  of  quasi-streamwise  vortices 
impinging  on  the  wall. 

Modelling  the  MS  pressure  depends  upon  the  vertical  velocity  held.  It  was 
shown  that  the  streamwise  spectra  of  v,  normalized  by  its  mean  square,  is  constant 
across  the  channel.  The  peak  of  the  spanwise  spectra,  however,  shifts  towards  lower 
wavenumbers  with  distance  from  the  wall.  The  location  of  the  peak  varies  as 
^0,  y-i/3  and  y~^  in  the  viscous  shear-layer,  buffer  layer,  and  logarithmic  region, 
respectively.  Further  efforts  are  necessary  to  scale  the  amplitude  of  the  spanwise 
V  spectra.  For  yly'  <  1.0,  the  broadband  correlations  were  shown  to  collapse  to  a 
single  curve  when  plotted  versus  y/y'. 
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Chapter  1 


Introduction 

1.1  Previous  Work 

Wall  pressure  fluctuations  (WPFs)  are  caused  by  the  passage  of  turbulent  flow  over 
the  wall  bounding  the  flow.  They  are  the  cause  of  flow-induced  structural  noise, 
which  can  cause  undesireable  airplane  cabin  noise  and  interference  with  acoustic 
sensors  and  are  also  important  for  wind-wave  generation  on  air-water  interfaces 
and  scouring  of  the  sea-bed  [Blake  (1986),  Willmarth  (1975)].  WPFs  are  related 
to  the  turbulent  flow  through  a  Poisson  equation,  and  they  are  a  function  of  all 
the  sources  in  the  flow  domain.  The  source  terms  consist  of  the  so-called  mean 
shear  (MS)  or  linear  term  and  turbulence-turbulence  (TT)  or  non-linear  terms. 
The  former  represents  the  interaction  between  the  mean  shear  and  the  fluctuating 
velocity  gradients,  while  the  latter  represent  the  interactions  between  fluctuating- 
velocity  gradients.  The  MS  term  is  also  called  the  “rapid”  term,  because  it  changes 
rapidly  in  response  to  the  flow  conditions,  whereas  the  TT  term  is  called  the  “slow” 
term  since  it  is  affected  only  through  non-linear  interactions  [Kim  (1989)].  Because 
of  the  complexity  of  the  three-dimensional  and  time-dependent  turbulent  flows,  it 
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has  been  difficult  to  determine  the  exact  relationship  between  velocity  held  sources 
and  WPFs.  The  present  state  of  knowledge  of  the  character  of  WPFs  and  their 
sources  is  due  to  the  contributions  of  many  researchers  over  the  past  40  years. 
Their  work  can  be  roughly  divided  into  three  areas:  modelling;  characterization 
in  frequency  and  frequency-wavenumber  space;  determination  of  source  terms  and 
their  relationships  to  WPFs,  in  physical  space.  In  this  review,  only  references  that 
directly  relate  to  the  current  work  will  be  mentioned.  Comprehensive  overviews  of 
the  work  prior  to  1988  can  be  found  in  Willmarth  (1975)  and  Eckelmann  (1988). 

Modelling 

Modelling  efforts  are  aimed  at  providing  a  WPF  model  in  wavenumber-frequency 
space  that  provides  the  forcing  function  for  flow-induced  structural  noise  predic¬ 
tions.  The  modal  response  of  a  structure  is  the  product  of  the  WPF  forcing 
function  and  the  structural  acceptances.  Corcos  (1963)  developed  a  model  that 
assumes  that  the  cross-spectral  density  of  pressure  is  the  product  of  three  func¬ 
tions;  the  autospectrum  of  the  point  pressure  fluctuations  and  two  exponential 
functions,  to  account  for  pressure  decorrelation  in  the  streamwise  and  spanwise 
directions.  The  Corcos  model  was  developed  to  describe  the  cross  spectrum  of 
pressures  dominated  by  the  convective  wavenumber  range.  As  initially  developed, 
it  has  questionable  validity  in  the  sub-convective  wavenumber  range  {kx  ^  u)/Uc, 
where  kx  is  the  streamwise  wavenumber,  C4  is  the  convection  velocity  and  oo  de¬ 
notes  the  radian  frequency)  because  measurements  and  applications  at  the  time 
were  from  the  convective  range.  In  many  subsequent  applications  [Blake  (1986)], 
the  underlying  wall  structure  has  its  major  acceptance  modes  in  the  subconvective 
wavenumber  range;  thus,  variants  of  the  Corcos  model,  which  attempt  to  model 
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the  subconvective  wavenumber  range  are  widely  used  for  flow-induced  structural 
noise  predictions  [Hwang  and  Maidanik  (1990),  Ko  and  Schloemer  (1992)]. 

Models  that  explicitly  include  the  source  terms  have  more  potential  for  accurate 
modelling.  Blake  (1971,  1986)  expressing  the  solution  to  Poisson  equation  in  a 
form  first  applied  by  Lilley  (1960)  obtained  a  double  integral  over  the  wall-normal 
coordinate,  which  involves  the  evaluation  of  the  five- dimensional  cross-spectral 
density  of  the  wall-normal  velocity  fluctuations,  <h22(|/,  |/^  K,  tc)  {y  is  the  wall- 
normal  coordinate,  K  denotes  the  planar  wavenumber  vector) .  He  further  assumed 
that  the  cross-spectral  density  can  be  separated  into  a  wall-normal  correlation 
function,  and  streamwise  and  spanwise  wavenumber  spectra  of  the  wall-normal 
velocity.  However,  even  these  much  simpler  functions  are  difficult  to  determine 
because  of  the  lack  of  data  that  may  be  used  to  define  the  spatial  structure  of  the 
wall-normal  velocity  field. 

Panton  and  Linebarger  (1974)  followed  in  this  vein  by  invoking  further  detail  in 
describing  the  variation  of  root-mean-square  (r.m.s.)  quantities  in  the  boundary- 
layer,  but  they  characterized  spatial  statistical  properties  of  the  sources  as  es¬ 
sentially  homogenous.  Recently,  Zawadski  et  al.  (1996)  evaluated  Blake’s  (1971) 
formulation  for  a  flow  with  adverse  pressure  gradient  and  separation  using  direct 
numerical  simulation  (DNS)  data.  They  were  able  to  define  some  of  the  essential 
correlations  and  predict  the  wall  pressure  spectra  which  had  the  correct  trends 
when  compared  with  DNS  pressure. 

Chase  (1980)  was  the  first  to  attempt  modeling  of  the  spatial  characteristics  of 
boundary-layer  structure  and  attempted  to  invoke  details  of  structure.  He  assumed 
that  the  wavenumber  spectra  of  the  wall-normal  velocity  scales  like  Ky  [K  denotes 
the  planar  wavenumber  amplitude),  based  on  the  conclusions  of  Morrison  and 
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Kronauer  (1969);  Chase  (1980)  then  obtained  a  form  for  both  the  MS  and  TT 
pressure  spectra  that  is  dependent  upon  eight  empirical  coefficients,  which  control 
the  shape  and  level  of  the  pressure  spectra. 

These  models  point  out  that  there  is  a  need  for  more  information  about  the 
source  terms,  and  more  information  about  the  wavenumber-frequency  character¬ 
istics  of  WPFs,  particularly  in  the  subconvective  wavenumber  range  range,  with 
which  to  validate  the  models. 

Characterization  in  spectral  space 

Analysis  of  the  frequency  spectra  of  WPFs  at  a  point  on  the  plane  was  perhaps 
the  first  tool  used  to  study  the  nature  of  the  source  terms.  By  associating  var¬ 
ious  frequency  ranges  with  regions  of  the  boundary-layer  it  has  been  possible  to 
hypothesize  where  in  the  boundary-layer  the  dominant  sources  are  located.  In  the 
ensuing  discussions  frequency  and  streamwise  wavenumber,  kx  (=  oj/Uc)  will  be 
used  interchangeably,  since  in  the  literature  both  appear,  and  can  be  related  to 
each  other  through  Taylor’s  frozen  convection  hypothesis.  Accordingly,  most  early 
work  in  source  modelling  addresses  principally  convected  pressure  fluctuations. 

Bradshaw  (1967)  argued  that  between  kxS  4  (5  is  the  boundary-layer  thick¬ 
ness)  and  a  wavenumber  inversely  proportional  to  the  height  of  the  viscous  shear 
layer,  the  spectra  should  have  a  decay  rate.  This  is  called  the  “overlap  re¬ 
gion”  because  the  pressure  spectra  will  collapse  for  varying  Reynolds  numbers 
when  normalized  both  by  outer  and  inner  variables.  It  has  been  shown  by  many 
researchers  [Farabee  and  Casarella  (1991),  Keith  et  al.  (1992),  Olivero-Bally  et 
al.  (1993)]  that  there  is  a  region  with  a  k~^  decay,  which  increases  in  size  with 
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Reynolds  number  [Panton  (1974),  Blake  (1986)].  Kraichnan  (1956),  considering 
just  the  Poisson  equation  solution  for  the  MS  pressure,  determined  that  in  the 
low  wavenumbers  the  spectra  should  have  a  k'^  behavior.  This  was  conhrmed  by 
Farabee  and  Casarella  (1991),  who  used  noise-cancellation  techniques  to  obtain 
accurate  low-frequency  spectra.  At  high  wavenumbers,  a  k~^  behavior  is  expected. 
If  inner-scaling  is  used,  the  k~^  behavior  appears  universal  (ie.,  independent  of 
Reynolds  number)  [Farabee  and  Casarella  (1991),  Keith  et  al.  (1992),  Olivero- 
Bally  et  al.  (1993),  Chang  et  al.  (1994)]. 

In  an  effort  to  shed  some  light  on  the  relationship  between  wavenumber  spec¬ 
tra  and  source  locations.  Panton  and  Linebarger  (1974),  starting  with  the  formu¬ 
lation  of  Kraichnan  (1956),  obtained  an  approximate  mimerical  solution  for  the 
contribution  of  the  MS  term  to  the  pressure.  They  validated  the  k^^  decay  rate 
and,  in  addition,  showed  the  influence  that  various  regions  of  the  boundary-layer 
have  on  the  streamwise-wavenumber  pressure  spectra.  They  demonstrated  that 
MS  sources  from  the  outer  region  (y  >  0.25)  affected  only  the  low-wavenumbers 
{kxS  <  20),  up  to  the  location  where  the  k~^  decay  rate  began.  The  middle 
{SS.du/ur  <  y  <  0.25)  {p  is  the  kinematic  viscosity,  and  Ur  is  the  friction  velocity) 
and  inner  regions  (y+  <  33.6)*  affected  the  higher-wavenumbers.  They  also  showed 
that  sources  from  the  middle  region  become  dominant  at  the  lowest  wavenumbers, 
hinting  at  the  possibility  that  all  regions  of  the  boundary-layer  contribute  to  the 
low- wavenumber  pressure. 

Blake  (1986)  showed  that  the  k~^  and  k~^  behaviors  can  be  obtained  by  just 
considering  the  MS  term,  and  making  estimates  for  the  mean  velocity  derivative 
and  the  cross-spectral  density  of  the  wall-normal  velocity  fluctuations  appropriate 

*y+  =  yUr/v 
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to  the  various  regions  of  the  boundary-layer  (viscous  shear-layer,  logarithmic  and 
wake  regions).  The  fact  that  experimental  data  conhrms  these  behaviors  lends 
credence  to  the  belief  that  the  MS  term  is  more  important  than  the  TT  terms. 
It  also  conhrms  the  conjecture  that  there  is  correlation  between  certain  regions  of 
the  boundary-layer  wavenumber  ranges.  In  summary,  the  weight  of  the  evidence 
is  that  the  k'^  range  is  due  to  the  wake  region;  the  k~^  is  due  to  the  logarithmic 
region;  the  k~^  range  is  due  to  the  viscous  shear-layer. 

Characterization  in  physical  space 

Perhaps  the  hrst  experiments  that  showed  the  two-dimensional  physical  charac¬ 
teristics  of  wall  pressure  were  performed  by  Emmerling,  Meier  and  Dinkelacker 
(1973);  using  optical  interferometry  techniques  they  were  able  to  view  the  flexi¬ 
ble  wall  deformations  caused  by  the  wall  pressure  fluctuations  beneath  a  turbu¬ 
lent  boundary  layer.  It  is  now  thought  that  coherent  turbulent  structures  are 
very  important  for  wall  pressure  generation  and  that  high-amplitude  wall-pressure 
events  contribute  a  significant  proportion  to  the  mean-square  wall-pressure,  even 
though  they  only  occur  a  fraction  of  the  time  [Karangelen  (1991),  Abraham  and 
Keith  (1995)].  Based  on  experimental  and  numerical  evidence,  it  has  long  been 
known  that  the  boundary- layer  is  populated  by  coherent  events.  Overviews  of 
the  key  developments  in  coherent-event  research  are  given  by  Cantwell  (1981)  and 
Robinson  (1991).  When  Cantwell  (1981)  wrote  his  review  article,  the  accepted 
boundary-layer  events  were:  (1)  low-speed  sublayer  streaks,  possibly  related  to 
pairs  of  streamwise  vortices;  (2)  the  lift-up  of  streaks  and  their  turbulent  burst¬ 
ing,  occuring  at  an  outer  scaled  frequency,  pointing  to  inner-layer  bursting  events 
controlled  by  outer  flow  events;  (3)  tilted  shear  layers,  the  upstream  sides  of  which 
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were  sites  of  intense  turbulence  activity;  (4)  the  inner  buffer  region  as  the  site  for 
most  of  the  turbulence  production.  Based  on  these  notions,  researchers  studied  the 
relationship  between  velocity  and  wall  pressure  events  using  conditional  sampling 
techniques.  Thomas  and  Bull  (1983)  used  VITA  techniques  in  an  experimental 
turbulent  boundary-layer  to  obtain  conditional  samples  of  pressure  traces  related 
to  high- amplitude  velocity  perturbations.  They  concluded  that  high- amplitude 
wall-pressure  events  are  produced  by  near-wall  shear  layers  and  that  such  events 
consist  of  a  region  of  positive  pressure  with  upstream  and  downstream  regions  of 
negative  pressure.  Johansson  et  al.  (1987),  also  using  conditional  sampling  tech¬ 
niques,  showed  that  near-wall  shear  layers  and  positive  pressure  events  have  a 
bi-directional  relationship  (i.e.,  the  same  events  are  detected  whether  velocity  or 
pressure  is  used  as  the  detection  criteria)  and  hypothesized  that  the  MS  term  was 
primarily  responsible  for  the  high-amplitude  wall-pressure  events. 

The  1980’s  and  90’s  saw  the  increasing  usage  of  numerical  simulations  of  tur¬ 
bulent  channel  and  boundary- layer  flows.  Simulations  allow  access  to  the  three- 
dimensional  and  temporal  velocity  held  data,  from  which  coherent  events  can  be 
extracted  more  easily.  Robinson  et  al.  (1988,  1989)  and  Robinson  (1991a,b) 
used  a  database  obtained  from  a  direct  numerical  simulation  (DNS)  of  a  low- 
Reynolds  number  boundary-layer,  determined  that  quasi-streamwise  vortices  and 
near-wall  shear  layers  are  among  the  primary  turbulent  structures  in  wall-bounded 
flows,  and  that  spanwise  vortices  in  the  logarithmic  region  occur  less  frequently. 
They  hypothesized  that  horseshoe  and  hairpin  vortices  proposed  by  many  re¬ 
searchers  [Cantwell  (1981),  Robinson (199 lb)  ],  which  are  a  combination  of  quasi- 
streamwise  vortices  and  spanwise  vortices,  were  less  common  than  the  individual 
vortices.  They  hypothesized  that  near-wall  shear  layers  are  closely  related  to  quasi- 
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streamwise  and  spanwise  vortices;  they  are  thought  to  form  when  high-speed  outer- 
region  fluid  collides  with  the  slower-moving  fluid  ejected  from  the  upward-directed 
sides  of  the  quasi-streamwise  or  spanwise  vortices  [Robinson  et  al.  (1988)]. 

Visualization  of  turbulence  structures  showed  that  quasi-streamwise  vortices 
generate  regions  of  negative  wall-pressure  beneath  low-pressure  vortex  cores  and 
ejections,  while  positive  pressures  have  been  shown  to  occur  beneath  sweeps  [Robin¬ 
son  et  al.  (1988),  Bernard  and  Thomas  (1993)].  Using  the  DNS  of  Kim  et 
al.  (1987),  Johansson,  Alfredsson,  and  Kim  (1988)  conditionally  averaged  the 
streamwise  velocity  to  obtain  a  three-dimensional  and  temporal  view  of  near-wall 
shear  layers  and  their  relationship  with  wall  pressure.  They  found  convincing 
evidence  that  high-amplitude  wall-pressure  events  are  related  to  near-wall  shear 
layers;  the  pressure  footprint  they  obtained  seems  consistent  with  the  findings  of 
Thomas  and  Bull  (1983).  Using  the  DNS  database  of  Choi  and  Moin  (1990), 
Lueptow,  (1995)  detected  pressure-gradient  events,  and  found  a  pressure  footprint 
consistent  with  previous  workers  [Thomas  and  Bull  (1983),  Johansson,  Alfreds¬ 
son,  and  Kim  (1988),  Wilczynski  (1992)].  Using  the  same  database,  Abraham  and 
Keith  (1995)  detected  high-amplitude  negative  and  positive  wall-pressure  events 
and  computed  their  planar  wavenumber  spectra. 

Casarella  and  his  coworkers  have  pursued  the  relationships  between  coher¬ 
ent  events  and  wall  pressure  in  a  low-noise  experimental  facility  described  in 
Farabee  (1986);  Wilczynski  (1992)  found  that  the  relationship  between  wall  pres¬ 
sure  gradients  and  Reynolds  stress  events  was  stronger  than  that  between  pressure 
peaks  and  velocity  peaks,  and  detected  wall-pressure  gradient  events  and  related 
them  to  u  —  v  velocity  fields.  Kammeyer  (1995)  and  Russell  (1997)  used  pressure 
itself  to  detect  high-amplitude  wall-pressure  events,  and  conditionally  averaged 


the  associated  velocity  fields.  Russell  (1997)  used  band-pass  frequency  filtering 
to  obtain  pressure- velocity  fields  in  frequency  bands  commonly  associated  with 
inner-,  mixed-  and  outer-scaling.  His  results  are  generally  in  agreement  with  the 
notion  that  the  shear  layers  generate  high-amplitude  wall-pressure  events.  Kam- 
meyer  (1995),  using  wavelet  transforms,  found  similar  structures. 

Another  advantage  of  numerical  simulations  is  that,  from  the  two-dimensional 
wall  pressure  fields,  the  frequency- wavenumber  (three-spectra)  can  be  obtained, 
although  their  computation  is  quite  expensive  [Handler  et  al.  (1984),  Chang  et 
al.  (1994)].  A  long  time  series  is  necessary  to  obtain  a  reasonably  well- resolved, 
yet  smooth,  spectrum;  grid  resolution  must  match  the  frequency  resolution  in  or¬ 
der  to  resolve  the  convective  ridge  and  obtain  the  correct  high  frequency  behavior. 
Choi  and  Moin  (1990)  computed  a  well-resolved  wavenumber-frequency  spectra 
using  a  simulation  of  fully  developed  channel  flow.  They  were  able  to  evaluate 
various  methods  of  computing  convection  velocity  and  study  the  similarity  form  of 
the  frequency  spectrum.  However,  there  were  questions  about  the  behavior  of  their 
spectra  in  the  subconvective  wavenumber  range.  Chang  et  al.  (1994)  compared  the 
wavenumber-frequency  spectra  computed  using  LES  and  DNS  and  studied  their 
characteristics,  particularly  in  the  sub-convective  wavenumber  range.  They  ob¬ 
tained  much  lower  subconvective  wavenumber  energy  than  Choi  and  Moin  (1990). 
Singer  (1996)  performed  an  LES  of  a  turbulent  boundary-layer,  and  obtained  the 
frequency  spectra  of  wall  pressure,  at  a  Reynolds  number  5  times  higher  than  that 
of  Chang  et  al.  (1994). 

Numerical  simulations  have  also  allowed,  for  the  first  time,  analysis  of  the  source 
terms.  Kim  and  Lee  (1989)  and  Kim  (1989)  computed  the  wall-normal  distribution 
of  the  mean-squares  of  the  MS  and  TT  source  terms.  Their  data  showed  that  the 
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dominant  source  term  is  the  one  most  strongly  associated  with  streamwise  vortices. 

Dhanak  et  al.  (1997)  numerically  simulated  single  and  pairs  of  streamwise 
vortices  and  studied  the  space-time  characteristics  of  the  related  wall  pressure. 
They  showed  that  the  time-space  trace  of  the  pressure  footprint  closely  resembled 
that  detected  experimentally.  They  obtained  frequency  and  spanwise  wavenumber 
power  spectra,  which  showed  the  same  roll-offs  obtained  for  the  total  pressure, 
although  their  spectra  levels  were  functions  of  an  externally  imposed  flow  field  and 
vortex  parameters. 

In  summary,  accurate  WPF  modelling  suffers  from  lack  of  information  about 
the  source  terms  and  the  well-resolved  wavenumber-frequency  pressure  spectra  that 
can  be  used  for  validation.  The  weight  of  the  research  into  turbulent  structures 
has  lead  to  the  conclusion  that,  at  least  in  low  Reynolds  number  flows,  quasi- 
streamwise  vortices  and  near-wall  shear  layers  are  the  dominant  coherent  events. 
Turbulence  simulation  databases  have  been  used  to  study  in  depth  the  nature  of 
wall  pressure,  but  have  been  used  only  sparingly  to  study  the  nature  of  the  source 
terms. 


1.2  Objectives 

The  objective  of  this  work  is  to  use  a  DNS  database  to  study  the  relationships 
between  the  velocity  held  source  terms  and  the  wall  pressure.  The  approach  used 
here  is  to  determine  which  source  terms  are  important  for  wall  pressure  generation, 
and  then  to  correlate  these  dominant  source  terms  to  physical  events  in  the  flow- 
field.  This  will  be  accomplished  by  identifying  the  regions  (buffer  layer,  logarithmic 
region  etc.)  that  give  the  most  significant  contributions  to  the  pressure  spectra 
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in  various  wavenumber  ranges,  and  then  by  relating  the  dominant  terms  to  the 
turbulent  structures  in  the  boundary-layer. 

Chapter  2  describes  the  mathematical  formulation  and  numerical  techniques  of 
the  turbulence  simulations;  it  contains  the  derivation  of  the  Poisson  equation  for 
wall  pressure  and  its  analytical  solution  for  a  channel  flow.  Chapter  3  validates 
the  velocity  and  wall  pressure  helds,  comparing  statistics,  spectra  and  correlations 
to  published  data.  In  Chapter  4  the  total  wall  pressure  is  decomposed  into  partial 
pressures  due  to  individual  terms  and  regions  and  the  dominant  partial  pressures 
are  determined.  Chapter  5  is  concerned  with  the  characteristics  of  the  sources 
associated  with  the  dominant  partial  pressures;  it  is  found  that  these  source  terms 
are  closely  related  to  coherent  turbulent  events.  Chapter  6  investigates  modelling 
of  the  MS  pressure  and  uses  the  channel  flow  data  to  begin  to  model  the  terms  in 
the  MS  pressure  formulation.  Finally,  concluding  remarks  and  recommendations 
for  further  research  can  be  found  in  Chapter  7. 
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Chapter  2 


Mathematical  Formulation 


This  chapter  describes,  in  detail,  the  mathematical  formulation  and  numerical 
techniques  used  in  the  turbulence  simulations  of  the  fully  developed  turbulent 
channel  flow  that  was  used  to  generate  the  pressure  databases.  Also  it  describes 
the  derivation  and  solution  of  the  Poisson  equation  for  fluctuating  wall  pressure. 
The  governing  equations  are  described  in  Section  2.1;  the  numerical  techniques  are 
presented  in  Section  2.2;  mathematical  dififerences  between  the  modified  pressure 
obtained  directly  from  the  simulations  and  the  physical  pressure  are  derived  in 
Section  2.4.  Derivation  of  the  Poisson  equation  for  fluctuating  pressure  is  described 
in  Section  2.5  and  its  analytical  solution  are  described  in  Section  2.6. 


2.1  Governing  equations 


The  time-dependent,  three-dimensional  Navier-Stokes  and  continuity  equations  for 
incompressible  flow  are 


dui  d  dP  1  d'^Ui 

dt  dxj  dxi  ^  Ref,  dxjdxj 


(2.1.1) 
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duj 

dxi 


0. 


(2.1.2) 


The  three  velocity  components  of  the  large-scale  motions  in  the  three  directions* 
are  u,  v  and  vj  or  ui,  U2  and  P  is  the  physical  pressure,  whereas  the  modihed 
pressure  (see  below)  is  denoted  by  (p.  The  flow  variables  are  made  dimensionless 
by  the  half-channel  width,  6,  a  reference  velocity  Ui  and  the  kinematic  viscosity, 
u.  The  reference  Reynolds  number  Rci  is  deflned  as 

Re^  =  ^.  (2.1.3) 


2.2  Numerical  technique 

The  time-dependent,  incompressible  Navier-Stokes  equations  (2.1.1)  in  vector  form 
can  be  written  as 

+  h(u)  = -VP  +  uW^u  (2.2.1) 

at 

where  h{u)  is  the  nonlinear  term.  The  nonlinear  term  h(u)  can  be  written  in  the 
skew-symmetric  form 

h{u)  =  -^  [u  •  Vu  +  V  •  (uu)]  (2.2.2) 

or  in  the  rotational  form 

/i(u)  =  — u  X  a;  +  •  (uu)  (2.2.3) 

where  a;  =  V  x  u  is  the  vorticity  vector.  The  term  | V  •  (uu)  goes  into  the  modified 
pressure  term  if  the  rotational  form  is  used,  and  the  resultant  computed  pressure 
is  P  +  |V  ■  (uu).  Zang  (1991)  showed  that  the  skew-symmetric  form  is  the  most 

*The  streamwise  direction  will  be  denoted  by  x  or  xi,  the  vertical  direction  by  y  or  X2,  and 
the  spanwise  direction  by  2:  or  0:3. 
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accurate,  while  the  rotational  form  is  less  expensive.  Both  forms  have  been  used 
in  the  present  simulations. 

The  governing  equations  have  been  advanced  in  time  by  a  fractional  time  step 
method  [Zang  and  Hussaini  (1986)].  Fourier  expansions  were  used  in  the  stream- 
wise  and  spanwise  directions,  with  a  Chebychev  collocation  method  in  the  wall- 
normal  direction. 


Temporal  discretization 


The  Navier-Stokes  eqirations  (2.2.1)  can  be  written  in  the  form 


du  d'^u 

__  -V 


where 


H  = -iu  .  Vu  -  iv  .  (uu)  +  ^  iv^u 


(2.2.4) 


(2.2.5) 


In  the  fractional  time  step  method,  an  intermediate  velocity  u*  is  computed  first 
from  a  Helmholtz  equation  (the  momentum  equations  without  the  pressure  term); 
then,  a  pressure  correction  is  obtained  by  computing  a  Poisson  equation  for  the 
modified  pressure  dehned  in  (2.4.10)  to  enforce  the  continuity  constraint.  Equation 
(2.2.4)  is  discretized  as 

*  n  r) 

(2.2.6) 


u*  -  u” 


At 


=  H  +  i/ 


dy^ 


and 


n+l  _  * 

,  =  -vr+h 

At  ^ 


(2.2.7) 

Taking  the  divergence  of  the  latter  equation  and  using  the  divergence- free  condition 
for  yields  a  Poisson  equation  for  the  modified  pressure 


At 


(2.2.8) 
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The  intermediate  velocity  u*  and  modified  pressure  0  are  computed  from  (2.2.6) 
and  (2.2.8);  subsequently,  (2.2.7)  gives  velocity 

The  vertical  diffusion  term  in  equation  (2.2.6)  is  treated  by  the  Crank- 
Nicholson  scheme,  which  is  implicit  and  second  order  accurate  in  time.  The  nonlin¬ 
ear  term  H,  which  includes  the  nonlinear  term  and  the  remaining  diffusion  terms, 
is  advanced  by  using  a  low  storage  third-order  Runge-Kutta  method  [Williamson, 
1980].  Thus,  the  time  advancement  process  of  the  governing  equations  is 


STAGE  1 


1  - 


u*'  = 


1  + 

|AtVV  =  V-u*' 
u'  = 


u”  +  |AtH" 


(2.2.9) 


STAGE  2 


1  — 

24  oy^ 


u*"  = 


+i|At(H'-  |H” 


^AtVV'  =  V-u 


u"  =  u*"-^AtV0" 


(2.2.10) 


STAGE  3 


1  - 


u 


*n+l  _ 


1 


dy^ 


u 


^rr 


_L  o  A  154  TTf  I  153  5  TTr 


iAtVV^^  =  v-u*”+^ 

u"+i  =  u*”+^  -  iAtV0”+^ 
these  equations  are  solved  in  Fourier  space. 


(2.2.11) 


The  CFL  value  for  the  RK3  scheme  is  computed  from 


CEL  =  max 


At 


u 


Ax 


+ 


Ay 


+ 


w 


Az 


(2.2.12) 
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For  numerical  stability,  the  RK3  scheme  requires  CFL  <  0.55.  The  value  used  in 
the  simulations  was  between  0.2  and  0.4. 

Spatial  discretization 

In  the  simulations  Fourier  expansions  are  used  in  the  streamwise  and  spanwise 
directions  and  Chebychev  expansions  are  employed  in  the  wall-normal  direction. 
The  Fourier  series  expansions  of  velocity  are 

fhL-l 

u(x,  y,  2:,  f)  =  ^  ^  (2.2.13) 

k  k  —IE. 

i\>x  —  2  —  2 

where  k^  =  ^k^,  k^  =  ^k^,  and  and  are  the  imposed  periodicity  lengths 
in  the  streamwise  and  spanwise  directions;  Nx  and  are  the  number  of  points 
in  the  streamwise  and  spanwise  directions.  The  collocation  points  in  the  periodic 
directions  are 


2Tri 

NxLx' 

i  =  0,l,.. 

• ,  -  1 

(2.2.14) 

2TTk 

k  =  0,1, . 

(2.2.15) 

The  Chebychev  expansion  of  the  Fourier  coefficient  u(A;i,y,  kz^t)  is 

Ny 

u{kx,y,kz,t)  =  ^u{kx,j,kz,t)Tj{^)  (2.2.16) 

j=0 

where  Ny  is  the  number  of  grid  points,  j  is  the  Chebychev  collocation  point,  and 
Tj{^)  is  the  Chebychev  polynomial, 

r,(0  =  cos(jcos“^0>  j  =  0,  l,...,iVj;.  (2.2.17) 

In  the  wall-normal  direction,  the  velocities  are  defined  at  the  full-grid  points 

Cy=cos^,  j  =  0,l,...,Ny  (2.2.18) 

I\y 
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and  the  pressure  is  defined  at  the  staggered  (half-)  grid  points 

e.+i  =  cos  j  =  0,l,...,Ny-l  (2.2.19) 

I\y 

The  momentum  equations  are  enforced  at  the  full-grid  points  and  the  continuity 
equation  is  enforced  at  the  half-grid  points. 

Matrix  inversion 

Since  the  Laplacian  operator  in  Fourier  space  takes  the  form 

=  -kl  -  kl  +  (2.2.20) 

the  left  sides  of  equations  (2.2.9-2.2.11)  are  linear  functions  of  the  second  derivative 
operator  Applying  the  Chebychev  expansions  (2.2.16)  to  this  operator  results 
in  a  linear  full  Ny  x  Ny  matrix,  and  the  resultant  forms  of  equations  (2.2.9-2.2.11) 
can  be  expressed  in  a  general  form 


LU  =  F.  (2.2.21) 

The  inversion  of  the  matrix  L  requires  a  large  amount  of  computing  time  and  ad¬ 
ditional  storage.  To  improve  the  computing  efficiency,  the  preconditioning  matrix 
method — an  iterative  technique  with  fast  convergence — is  used  in  the  simulation. 

If  U°  is  an  initial  guess  of  U,  the  initial  residual  is 

R°  =  F  -  LU°  (2.2.22) 

which  can  be  also  written  as 

L  (U-U°)  =  R°.  (2.2.23) 
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Applying  a  finite  difference  discretization  technique  to  equations  (2.2.9-2.2.11)  re¬ 
sults  in  a  finite  difference  preconditioning  matrix  G, 

G  AU°  =  R°  (2.2.24) 


in  which  AU°  is  the  finite  difference  solution  of  U  —  U*^.  The  iterative  direc¬ 
tion  is  adjusted  by  the  Conjugate  Residual  method  proposed  by  Hestenes  and 
Stiefel  (1952)  to  improve  the  convergence  rates: 

Ufc+i  =  +  cifcAU^  (2.2.25) 


with 

(L  AU^,  i?^) 

“  (laumTmj^' 

Then,  the  residual  and  the  iteration  equation  at  the  fc-|-l-th  step  are 

R^+i  =  R^  -  ttfcL  AU^ 


(2.2.26) 


(2.2.27) 


and 

GAU^+i  =  R^+b  (2.2.28) 

This  process  continues  until  convergence  is  reached.  The  preconditioning  matrix  G 

used  above  is  obtained  by  performing  a  second-order  finite  difference  approximation 

to  the  operator  L  at  the  Chebychev  collocation  points 

—2  2  —2 

hj-i{hj  +  hj{hj  -h 

j  =  l,...,iV',-l  (2.2.29) 

uo  =  0,  (2.2.30) 

UN,  =  0  (2.2.31) 

where  hj  =  yj  —  yj+i-  This  process  shows  an  extremely  fast  convergence  rate 
and  only  takes  less  than  30  iteration  steps  to  reach  convergence  in  the  present 
simulations. 
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Dealiasing 


The  nonlinear  terms  of  the  Navier-Stokes  equations  and  the  velocity  held  source 
terms,  presented  later  are  subject  to  aliasing  errors  if  measures  are  not  taken  to 
eliminate  them.  The  2/3  dealiasing  rule  was  used  to  eliminate  aliasing  errors  in  the 
computation  of  the  nonlinear  products.  Aliasing  occurs  because  the  multiplication 
of  two  functions  with  M  resolved  components  results  in  a  function  with  energy  in 
2M  components,  with  the  highest  M  modes  unresolved.  Fourier  transforming  the 
product  will  result  in  the  unresolved  energy  being  aliased  back  into  the  resolved 
modes.  In  the  pseudo-spectral  method  used  here,  the  nonlinear  terms  are  computed 
in  physical  space,  but  need  to  be  transformed  into  Fourier  space  for  solution.  In 
the  2/3  dealiasing  rule  the  entire  simulation  is  run  on  a  grid  which  has  3Mx/2 
and  3M2/2  collocation  points  in  each  of  the  planar  directions,  respectively.  The 
following  dealiasing  method  is  described  in  terms  of  x,  but  the  concepts  are  the 
same  for  .2  and  x  —  z.  Starting  with  velocity  fields  with  only  resolved  modes, 
multiplication  will  result  in  a  function  with  energy  in  2Mx  modes.  However,  since 
the  solution  is  on  a  grid  with  modes,  the  aliased  energy  will  be  folded  about 
the  3Mx/2  +  1  mode  so  that  the  aliased  energy  will  end  up  only  in  the  top  1/3 
modes,  M^.  -|-  1  <  |ml  <  3Mx/2.  The  top  1/3  modes,  which  would  contain  the 
aliasing  errors  are  set  to  zero.  The  resulting  velocity  fields  have  dealiased 
modes,  or  2/3  of  the  modes  upon  which  the  solution  is  computed. 

Boundary  conditions 

No-slip  and  no-flow-through  boundary  conditions  are  enforced  at  the  walls.  Peri¬ 
odic  boundary  conditions  are  enforced  at  the  streamwise  and  spanwise  boundaries. 
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Since  the  method  of  the  fractional  time  step  introduces  the  intermediate  velocity 
u*,  its  boundary  conditions  also  need  to  be  defined.  It  can  be  expressed  as 

u*  =  (2.2.32) 

from  (2.2.7).  is  approximated  by  the  second  order  accurate  Taylor  expansion 
at  the  boundaries  and  yields  a  third  order  accurate  approximation  for  u*; 

u*  =  +  OiAt^)  (2.2.33) 

which  has  the  same  order  accuracy  as  the  RK3  scheme.  The  boundary  condition 
for  u*  is  written  as 

s 

u*  =  AtS/cfA  +  At2_V</)”  (2.2.34) 

(-/  T 

at  the  walls. 


2.3  Large-eddy  simulation 


In  addition  to  the  DNS  that  produced  most  of  the  data  presented  in  this  work,  some 
large-eddy  simulations  (LES)  were  also  run  and  will  be  discussed  in  Chapter  3.  In 
large-eddy  simulations  the  flow  variables  are  decomposed  into  a  large-scale  (or 
resolved)  component,  denoted  by  an  overbar,  and  a  subgrid-scale  component.  The 
large-scale  component  is  defined  by  the  Altering  operation. 


f{Xi,X2,Xs,t)  =  )f{x[,  X2,x'^,  t)dx[dx2dx/^ 


Id 


(2.3.35) 


i=l 


where  D  is  the  computation  domain  and  Gi  is  the  filter  function  in  the  i-th  direc¬ 
tion.  Applying  the  filtering  operation  to  the  Navier-Stokes  and  continuity  equa¬ 
tions  yields  the  governing  equations  for  the  large  scale  motions. 


dui  d 


dcf)  dTij  1  d'^u 


+ 


dxi  dxj  Ree  dxjdxj 


(2.3.36) 
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The  subgrid-scale  stress  term, 


Tij  'lljtlj 


(2.3.38) 


represents  the  contribution  from  the  subgrid  scales  and  must  be  modeled.  In 
this,  the  dynamic  model  [Germano  et  a/.(1991),  Lilly  (1992)]  has  been  applied. 
The  anisotropic  part  of  the  SGS  stresses  are  parametrized  by  an  eddy  viscosity 
assumption 


where  Sij  is  the  large-scale  strain  rate  tensor  is  given  by 

c..  =,  1 

2  \dxj  dxi )  ’ 

and  the  eddy  viscosity  ut  is  given  by 


(2.3.39) 


(2.3.40) 


ut  =  CAVI, 


(2.3.41) 


where  l^l  =  \j2SijSij  and  A  =  (dA^jA^^A^)^/^  is  the  filter  width  at  the  grid  scale. 
The  coefficient  C  is  determined  from  the  energy  content  of  the  smallest  resolved 
scales  using  the  least  squares  approach  of  Lilly  (1992), 


1  CjjMjj 


(2.3.42) 


Here,  Cij  =  UiUj  —  UiUj  “resolve  the  turbulent  stresses”  and  Mij  =  —  (5ij  where 

(3ij  =  A^|S'|S'j^-  and  aij  =  A  Filtering  at  a  test  scale  with  filter  width 

A  >  A  is  denoted  by  A 


The  numerical  solution  of  the  filtered  equations  is  carried  out  using  the  scheme 
described  before;  however,  the  viscosity  u  is  replaced  by  <  ut  >  (where  <  ur  > 
is  the  plane-averaged  turbulent  eddy  viscosity  computed  using  the  dynamic  model); 
the  remainder  of  the  SGS  stress  is  advanced  explicitly  using  the  RK3  method. 


2.4  Modified  versus  physical  pressure 


As  shown  in  Section  2.2,  the  modified  pressure,  0,  rather  than  the  physical  pressure, 
P,  is  solved  in  the  simulations.  In  this  section,  the  difference  between  the  two 
pressures  will  be  derived.  The  fractional  step  method  with  Crank-Nicholson  time 


advancement  yields  a  Helmholtz  equation 


=  -Mu")  +  ;t 


u  f  d'^u*  d'^u^' 


2  \  d'lf  dif 


The  pressure  correction  step  is  given  by 


(2.4.1) 


=  -V(jf 


(2.4.2) 


Taking  the  divergence  of  (2.4.2), 


V  •  -  V  •  u*  =  -AfV  •  (V</)" 


(2.4.3) 


and  observing  that  the  first  term  on  the  left-hand-side  is  zero  by  continuity,  gives 


a  Poisson  equation  for  (j), 


•  u* 

^  At 


(2.4.4) 


The  solution,  is  obtained  by  solving  (2.4.1)  for  u*,  then  (2.4.4)  for  and 
finally  (2.4.2). 


The  question  remains,  how  is  0  related  to  the  actual  pressure?  Add  the 


Helmholtz  and  correction  step  equations  (2.4.1)  and  (2.4.2)  together: 


u"+i  -  u" 


=  -fi(u”)  + 


u  /  d'^u” 

2  i  dy'^  ^  dy"^ 


(2.4.5) 


This  is  almost  the  momentum  equations  with  Crank-Nicholson  time  advancement 
except  for  that  fact  that  the  first  viscous  term  involves  u* .  The  viscous  term  with 
can  be  recovered  by  taking  the  second  derivative  in  y  of  (2.4.2) 


dy2  Qy2  +  ) 


(2.4.6) 
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Substitute  (2.4.6)  into  (2.4.5)  to  obtain 


- =  -h  +  - - 1 - 

At  ^  2[dy^  dy^ 


n+1 


+ 


pAt 


(2.4.7) 


dP 

dx. 


where  the  last  two  terms  must  equal  the  physical  pressure  term.  This  gives  the 
relationship  between  P  and  (j)  introduced  above  (if  the  rotational  form  is  used,  P 
still  contains  also  UiUi/2), 


=  X76-+1  _  ^ 

dxi  2  \  dy‘^  ) 


(2.4.8) 


Integration  of  the  three  equations  in  (2.4.8)  gives 

pn  +  l  _  _ 


-+1  _  M 


(2.4.9) 


2 

The  integration  yields,  at  most,  an  additive  constant  which  has  been  ignored  since 
we  are  interested  in  the  fluctuating  pressure,  and  the  constant  would  just  add  to 
the  mean  pressure.  When  using  the  rotational  form  of  the  nonlinear  terms,  the 
resolved  kinetic  energy  is  added  to  the  pressure  term.  Furthermore,  for 

subgrid-scale  models  which  only  models  the  anisotropic  part  of  the  subgrid-scale 
stress  tensor,  the  isotropic  part  ql  =  |  {txx  +  Tyy  -\-  Tzz)  is  added  to  the  pressure 


term. 


We  are  using  a  three  step  Runge-Kutta  procedure  for  explicit  time  advancement 
and  are  interested  in  obtaining  the  pressure  only  after  a  complete  time  step.  The 
effective  At  should  be  based  on  the  last  of  the  RK  time  increments.  The  third  RK 
step  advances  the  solution  from  4-  ^  to  +  At,  thus  the  time  increment  to 
use  in  (2.4.9)  is  yielding  the  pressure  equation 

+  i,.  =  •  (u"+‘u“+‘)  .  (2.4.10) 
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2.5  Derivation  of  Poisson  equation  for  fluctuat¬ 


ing  pressure 


The  Poisson  equation  for  fluctuating  pressure  is  obtained  by  taking  the  divergence 
of  the  incompressible  momentum  equations,  applying  Reynolds  decomposition  to 
separate  the  mean  and  fluctuating  quantities  and  then  applying  continuity.  If 
we  take  the  divergence  of  (2.1.1),  and  eliminate  the  viscous  terms  by  applying 
continuity  and  the  vector  identity. 


V  X  (V  X  U)  =  0, 


(2.5.1) 


the  Poisson  equation  for  the  total  pressure  results: 


d'^p 


=  -P 


d  ( duiUj 


(2.5.2) 


dxidxi  dxi  \  dxj  J 

We  now  decompose  the  total  pressure  and  velocities  into  their  mean  and  fluctuating 
quantities,  where  mean  quantities  are  denoted  by  capital  or  overbar,  the  fluctuating 
quantities  by  primes,  and  the  mean  is  a  temporal-  and  planar-average;  we  have 


p  =  P  +  p' 

u  =  U  +  vf  (2.5.3) 


The  spatial  and  temporal  averaging  satisfy  the  following  properties  characteristic 
of  Reynolds  averaging  operators: 


u'  =  0 

Uiu'j  =  0 

dui  dui 
dxj  dxj 


(2.5.4) 
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Expanding  (2.5.2)  we  obtain 


d^p  d^p' 
dxidxi  dxidxi 


d 


-P 


dxidxj 


(jJiUj  +  UiUj  +  UjU j  +  u^Uj^  . 


(2.5.5) 


We  take  the  mean  of  (2.5.5),  and  use  the  properties  (2.5.4)  to  obtain  a  Poisson 
equation  for  the  mean  pressure; 


d‘^P 

dxidxi 


Q  _ 


dxjdx 


(2.5.6) 


Now  substitute  (2.5.6)  into  (2.5.5)  to  obtain  the  Poisson  equation  for  the  fluctuat¬ 
ing  pressure, 

dxidxi  \  dxj  dxi 
where  the  first  term  on  the  right  hand  side  is  the  mean-shear  (MS)  (“linear”  or 
“rapid”)  source  term,  and  the  turbulence-turbulence  (TT)  (“nonlinear”  or  “slow”) 
source  term  is  given  by 


(2,5,7) 


rpTT 


dxidxj 


(«'«'  -  «'«')  . 


(2.5.8) 


The  TT  term  can  be  further  simplified  by  expanding  the  derivatives  of  the  terms 
within  the  parentheses  and  applying  continuity  to  obtain 


rpTT  _ 

dXj  dxi 


dxidXj 


u[u'y 


(2.5.9) 


Since  the  mean  velocity  for  a  channel  flow  is  Ui  =  ([/,  0,  0),  the  final  form  for  the 
Poisson  equation  for  fluctuating  pressure  in  a  channel  flow  is 


^^P  _  _  r rj^MS  I  rpTT 

dxidxi  ~  >- 


(2.5.10) 


The  MS  term  is  defined  as 


rj^MS 


dy  dx^ 


(2.5.11) 


25 


the  total  TT  term  is  the  summation  of  all  the  TT  terms, 

(2.5.12) 

i=l  j=l 

and  the  total  source  term  is 


T 


tot 


rpMS  pTT 


(2.5.13) 


For  a  plane  channel  the  in-plane  derivatives  of  the  planar  mean  terms  equal  zero 
and  the  TT  source  terms  become,  explicitly. 


ttt  _  f  du\ 

) 


TTT  _  rpTT  _  du  dv 
^12  -^21 

^tt  ^  ~ 


ij  dy^ 


TTT  _  rpTT  _  du  dw 
13  31  '^~3x. 


TTT  _  rpTT  _  dv  dw 
23  32  ^  '3^ 


T'rT  _ 
-‘-33  ~ 


dw 

~3z 


2 


(2.5.14) 


Note  that  the  primes  on  the  fluctuating  quantities  have  been  dropped  for  conve¬ 
nience. 


2.6  Solution  for  fluctuating  pressure 


In  this  section  we  will  find  a  solution  for  the  Poisson  equation  in  terms  of  the 
unknown  source  terms,  in  order  to  find  the  attenuation  of  the  source  terms  with 
wall-normal  distance  and  wavenumber,  and  the  influence  that  the  source  terms 
have  on  the  wall  pressure.  The  problem  that  we  would  like  to  solve  is 


d'^p 

dxidxi 


-T 


tot 


(2.6.1) 


where  . 
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Boundary  conditions 


Boundary  conditions  at  the  top  and  bottom  channel  walls  can  be  found  from  the 
wall- normal  Navier-Stokes  equation  for  the  total  velocities: 


Dx)  dp  1 

'm  ~  ^  ~ 


d'^v  d‘^v  d‘^v ' 

Re  \  dx^  dip  dz^ 


(2.6.2) 


The  LHS  of  the  equation  goes  to  zero  at  the  walls  because  of  the  no-slip  boundary 
conditions.  Since  v  =  0  everywhere  on  the  walls, 


d'^v 


This  leaves  the  condition 


dx'^ 


dp 

dy 


d^v 


y=±l 


dz^ 


=  0. 


(2.6.3) 


y=±l 


1  d'^v 


y=±l 


Re  dy"^ 


(2.6.4) 


y=±l 


Solution  method 

We  can  transform  the  Poisson  equation  for  the  fluctuating  pressure  into  Fourier 
space, 

where  the  denotes  a  planar  Fourier  transform;  K  =  +  k'l]  kx  =  ^  and 

kz  =  (2.6.5)  is  a  nonhomogenous,  linear,  second  order,  ordinary  differential 

equation,  which  can  be  solved  using  variation  of  parameters.  Note  that  there  is  no 
direct  time  dependence,  meaning  that  disturbances  in  the  flow  given  by  T  are  felt 
instantaneously  at  all  points  in  the  domain.  The  equation  is  elliptic;  the  pressure 
at  a  point  is  a  function  of  all  source  points  in  the  domain. 

The  complementary  solution  is  given  by 

Pc{y)  =  Aipi{y)  +  A2P2{y)  (2.6.6) 


27 


where  pi{y)  =  and  p2{y)  =  e  The  particular  solution  is  given  by 


hiy)  =  vi{y)h{y)  +  V2{y)h{y)- 

We  solve  for  the  particular  solutions  with 

fy  f*'°*'{r])p2{ri)dr] 


(2.6.7) 


vi{y)  =  -  j 


and 


/y 

- 

a. 


ao{y)W  [pi{y),p2{ri)\ 


(2.6.8) 


(2.6.9) 


io{y)W  [pi{y),P2{y)\ 

W  is  the  Wronskian  and  Uq  is  the  coefficient  of  the  second  order  term  in  (2.6.5). 
The  Wronskian  is  IT  =  —2K.  The  solution  becomes 


oKy_ 


Ky 


(2.6.10) 


p{y)  =  A^e^y+A2e-^y+\^-;^  j"  dp 

This  can  be  simplified  to  the  form 

p{y)  =  Aie^y  +  A2e~^y  +  ^  j  sinh  [K{y  -  p)]  dp.  (2.6.11) 


Application  of  boundary  conditions 


We  now  apply  the  boundary  conditions  (2.6.4).  We  can  set  the  lower  limit  of 
integration  to  anything  that  is  convenient,  as  a  long  as  we  apply  it  consistently. 
We  can  see  that  setting  it  to  ?/  =  —  1  (the  near- wall)  will  get  rid  of  the  integral  term 
when  we  apply  the  boundary  condition  at  y  =  —1.  Applying  the  near- wall  and 
far-wall  {y  =  1)  boundary  conditions,  respectively,  we  obtain  the  set  of  equations; 

AiKe-^  +  A2Ke^  =  G{-1) 


(2.6.12) 


AiKe^  +  A2Ke-^  =  G(l) 
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where 


and 


G(l)  = 


dp 

dy 


y=-l 

-Xi 


y=l 


Xi  =  4  sinh  [Kil  -  p)]  dp 

K  J-i 


(2.6.13) 


(2.6.14) 


Solution 


Finally,  it  can  be  shown  that  the  fluctuating  wall  pressure  can  be  expressed  as  two 
boundary  terms  and  an  integral  over  the  wall-normal  domain: 


p{K,  y  =  -l,p)  =  B+  -  B-  -  f 


where  G  is  the  Green’s  function, 


G{K,y=-Gp)  = 


^tot 


(K.,)G(A'.y  =  -l,r,)c(„,  (2,6,15) 


cosh[iF(l  —  p)] 


iFsinh(2/F)  ' 

The  near-  and  far-wall  upper  boundary  terms  are  given  by 


(2.6.16) 


B_=p_{K)i/{K,y  =  -l)  (2.6.17) 

B+^P4K)ff{K,y  =  l), 


respectively,  with  the  wavenumber-dependent  functions  given  by 


P-{K) 

MK) 


1 

iFtanh(2iF) 

1 

/Fsinh(2iF) 


(2.6.18) 
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(a)  y/s  (b) 

Figure  2.1:  Wall-normal  profiles  of  Green’s  function  for  wavenumber  {kx^kz)  com¬ 
binations:  - :  (1/3, 3/2); .  :  (1,3/2); .  :  (2,9/2); - :  (5,10); 

- :  (10.5,46.5)  [cutoff  for  LES]; - :  (31.67,70.5)  [cutoff  for  DNS],  (a)  Entire 

channel,  (b)  near-wall  region. 

Properties  of  the  Green’s  function  solution 

The  Green’s  function  for  fixed  y,  varies  with  rj  and  K.  Wall-normal  profiles  of  G 
for  K  in  the  range  occuring  in  the  DNS,  are  shown  in  Figure  2.1(a)  for  the  entire 
channel  and  in  the  near-wall  region  in  Figure  2.1(b).  Figure  2.1(a)  shows  that 
the  attenuation  in  the  lowest  wavenumber  is  only  about  12  dB  across  the  channel, 
meaning  that  the  lowest  wavenumber  sources  may  be  felt  all  the  way  across  the 
channel.  However,  the  attenuation  with  increasing  wavenumber  is  very  rapid,  as 
low  wavenumbers  e.y,  {kx  =  2,  =  9/2)  are  attenuated  by  47  dB  across  the 

channel.  Figure  2.1(b)  shows  the  profiles  out  into  the  logarithmic  region.  If  we 
classify  disturbances  as  being  insignificant  if  they  get  attenuated  by  20  dB  (i.e., 
they  become  1/100  of  their  original  value  at  y  =  y),  then  we  can  see  that  in  the 
buffer  layer  (say,  y+  =  18)  only  sources  with  wavenumbers  less  than  [kx  =  b,kx  = 
10)  are  significant. 
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Figure  2.2:  Wavenumber  profiles  of  boundary  terms:  —  :  /5_; . :  /5_  (high-K 

limiting  form); .  :  /?+; - '■  (^+  (high-K  limiting  form). 

Figure  2.2  shows  the  wavenumber-dependent  functions,  /5_(iF)  and  /5+(iF).  The 
effect  of  the  near-wall  boundary  term  on  the  pressure  at  the  near-wall  is  actually 
larger  than  0  dB  for  iT  <  1  and  then  falls  off  slowly  for  larger  values  of  K.  The 
far-wall  boundary  term  is  greater  than  0  dB  for  K  <  0.65,  and  then  decays  rapidly 
for  higher- KT.  The  high- wavenumber  limiting  forms  of  P{K)  are  given  by 

0-{K)  ~  ^  (2.6,19) 
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For  large  values  of  K,  G  becomes  a  simple  function; 

G{K,y  =  -l,rj)  =  (2.6.21) 

It  can  be  seen  that  G  is  singular  as  iF  — 0,  which  means  that  it  heavily  weights 
the  low-iF  boundary  layer  sources,  and  attenuates  the  high-iF  sources.  For  K  =  0 
Kim  (1989)  gives  the  formula 

G{0,y,y)  =  0.5iy-y).  (2.6.22) 
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Chapter  3 


Characteristics  of  Plane  Channel  Flow 


The  purpose  of  this  chapter  is  to  describe  the  numerical  simulations  and  to  show 
comparisons  between  pressure  and  velocity  data  with  those  from  the  literature. 
The  Chapter  is  divided  into  two  sections:  the  first  deals  with  validation  of  the 
velocity  fields  from  which  the  wall  pressure  and  velocity  field  source  terms  were 
computed;  the  second  deals  with  the  wall  pressure  fields  from  multiple  large  eddy 
and  direct  numerical  simulations. 


3.1  Velocity  fields 

In  the  course  of  this  work  multiple  simulations  were  run,  mainly  for  the  purpose 
of  studying  the  computational  aspects  of  computing  wall  pressure  wavenumber- 
frequency  spectra.  They  will  be  discussed  in  Section  3.2.  In  this  Section,  only  the 
velocity  fields  from  LES7  and  DNSS,  which  were  used  to  compute  the  partial  wall 
pressures  and  source  terms  in  Chapters  4  and  5,  respectively  and  for  modelling  of 
vertical  velocity  in  Chapter  6  will  be  discussed. 

In  these  simulations,  the  streamwise  domain  was  GttS,  50  percent  larger  than 
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LES7 

DNSS 

Kim  et  al  (1987) 

RCr  =  UtA/v 

171.8 

179.8 

~  180 

Ur/Uo  X  10^ 

5.265 

5.525 

5.49 

U})!  ti'f 

16.29 

15.57 

15.63 

UolUt, 

1.17 

1.16 

1.16 

5*/5 

0.1424 

0.1396 

0.141 

ei5 

0.0858 

0.0858 

0.087 

Cf  X  10'^ 

7.54 

8.25 

8.18 

Table  I;  Flow  parameters  for  turbulence  simulations 
that  used  by  Choi  and  Moin  (1990)  and  Kim  et  al.  (1987).  The  reason  for  the  larger 
domain  was  to  obtain  better  streamwise  wavenumber  resolution;  the  resolution  for 
LES7  and  DNSS  is  =  1/3,  (for  a  4n5  box  it  would  be  1/2);  the  maximum 
value  is  k^S  =  31.667  for  DNSS  and  11.667  for  LES7.  The  spanwise  domain  for  both 
simulations  is  4/37r(5.  The  spanwise  resolution  for  both  simulations  is  Ak^S  =  ?>/2^ 
with  maximum  values  of  k^S  =  70. S  for  DNSS  and  46.5  for  LES7. 

DNSS  had  192  x  96  x  96  collocation  points  in  the  streamwise,  spanwise  and 
wall-normal  directions,  respectively,  whereas  LES7  had  72  x  64  x  64  points  in  each 
of  the  directions.  The  streamwise  grid  resolution  for  DNSS  is  Ax^  =  17.6,  whereas 
it  is  44.8  for  LES7.  The  spanwise  grid  resolution  for  DNSS  is  Az^  =  7.8;  11.2  for 
LES7.  Piomelli  (1993)  showed  that  spanwise  resolution  Az'^  =  13  was  sufficient 
to  resolve  the  near- wall  streaky  structures  in  LES.  The  grid  point  nearest  the  wall 
was  Ay^  =  0.09  for  DNSS  and  0.21  for  the  LES7. 

Both  simulations  used  the  rotational  scheme  and  were  dealiased  in  the  stream- 
wise  and  spanwise  directions  using  the  2/3  dealiasing  rule  described  in  Section  2.2. 
Chang  et  al.  (1994)  showed  that  aliasing  errors  can  significantly  affect  the  high 
frequency  and  wavenumber  components  of  the  wall  pressure  spectra. 

Some  of  the  pertinent  simulation  parameters  are  listed  in  Table  I.  It  shows 
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(b)  y" 


Figure  3.1;  Mean  velocity  profiles.  . .  DNS5;  A  LES7.  (a)  Global  coordi¬ 

nates;  (b)  wall  coordinates. 

that  DNSS  and  Kim  et  al  (1987)  are  in  very  good  agreement.  LES7  has  a  lower 
value  of  Ur  which  is  attributable  to  decreased  grid  resolution  [Piomelli  (1993)]. 
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Mean  velocity  profiles 


Figure  3.1  shows  the  mean  streamwise  velocity  profiles  in  global  and  wall  coordi¬ 
nates.  In  Figure  3.1(a),  the  velocity  has  been  normalized  by  the  bulk  velocity, 

U,  =  ]^f_Udy,  (3.1.1) 

the  comparison  of  velocities  between  DNSS  and  LES7  is  very  good.  The  values 
of  Ub/ur  are  15.6  and  16.3  for  DNSS  and  LES7,  respectively,  in  good  comparison 
to  the  value  of  15.6  obtained  by  Kim  et  al.  (1987).  The  mean  velocity  profiles 
normalized  by  wall  coordinates  and  plotted  against  the  law  of  the  wall  ^ 

which  defines  the  viscous  shear-layer,  and  wT  =  2.5  In  y"*"  -|-  i7,  the  logarithmic 
region,  is  shown  in  Figure  3.1(b).  As  can  be  seen,  both  DNSS  and  LES7  agree 
very  well  in  the  viscous  shear-layer  but  diverge  slightly  in  the  buffer  layer  and 
logarithmic  region.  The  constant,  B,  is  5.5  for  DNSS  and  about  6  for  LES7.  The 
difference  in  intercepts  is  due  primarily  to  when  U  from  LES7  is  normalized 
by  Ur  from  DNS5,  it  collapses  to  the  DNSS  curve.  The  logarithmic  region  extends 
out  to  the  middle  of  the  channel  (y+  =  180)  with  no  trace  of  a  wake  region,  as 
expected  at  this  Reynolds  number. 

Two-point  correlation  coefficients 

Two-point  spatial  correlations  of  velocities  indicate  the  length  scales  of  the  turbu¬ 
lence  structures  within  the  channel.  Since  the  channel  flow  is  statistically  station¬ 
ary  and  homogenous  m.  x  —  z  planes,  the  spatial  correlations  are  only  dependent 
upon  the  spatial  separation  between  two  points  in  that  plane,  and  are  shown  only 
as  functions  of  two  variables;  a  streamwise  or  spanwise  separation,  (^  and  rj,  re- 


36 


Figure  3.2:  Two-point  velocity  correlations  with  streamwise  separations; -  Ruu] 

A  Ryy]  o  R^^.  (a)  ?/+  =  12,  (b)  =  150. 

spectively)  and  the  wall-normal  location,  y.  The  autocorrelations  are  defined  as 

<  u{x\  y)u{x  +  C;  y)  > 


Ruui^j  y)  — 


(3.1.2) 


<  yW  > 

Shown  in  Figure  3.2  are  the  streamwise  autocorrelations  obtained  from  LES7  at 
=  12  and  150.  They  show  that  u  has  the  largest  streamwise  correlations,  which 
is  probably  due  to  the  alternating  low-  and  high-  speed  streaks  in  the  near-wall 
region.  The  u  correlations  go  to  zero  with  a  separation  distance  of  about  =  6. 
This  indicates  that  the  the  longest  streamwise  structures  are  less  than  half  the 
box  length  {L^/S  =  Ott  ~  19)  and  that  the  streamwise  box-size  is  sufficient.  The 


37 


Figure  3.3;  Two-point  velocity  correlations  with  spanwise  separations; -  Ruu] 

A  o  (a)  =  12,  (b)  y+  =  150. 

spanwise  autocorrelations  of  the  velocities,  shown  in  Figure  3.3,  indicate  that  the 
spanwise  correlation  lengths  are  also  much  smaller  than  the  streamwise  lengths. 
The  Ruu  has  a  negative  peak  at  =  60  {rjl5  =  0.35)  which  could  be  indicative  of 
alternating  streaks.  These  compare  very  well  with  Kim  et  al.  (1987). 


The  streamwise  and  spanwise  correlations  of  wall  pressure,  are  shown  in  Fig¬ 
ure  3.4.  The  streamwise  correlations  go  to  zero  after  about  ^/S  =  2  —  3,  whereas 
the  spanwise  correlations  go  to  zero  very  slowly  and  do  not  reach  zero  for  the 


Figure  3.4:  Two-point  wall  pressure  correlations  with  (a)  streamwise  and  (b)  span- 
wise  separations. 

longest  spanwise  separations.  The  low  correlation  at  rj/S  =  Lz/2  (less  than  0.05) 
is,  however,  acceptable,  and  the  box  was  considered  sufficiently  wide. 
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Figure  3.5;  Shear  stress  profiles.  -  DNS5;  A  LES7 

where  S  and  u  have  been  used  for  normalization  and  <Ti2>  is  the  subgrid 
scale  stress,  defined  in  (2.3.38).  Thus,  for  a  steady  state  solution,  the  average 
total  shear  stress  converges  to  a  linear  function.  Figure  3.5  shows  that  the  total 
shear  stress  for  DNS5  and  LEST  are  linear,  indicating  that  the  statistical  samples 
are  well  converged.  Also  shown  in  Figure  3.5  are  the  Reynolds  stress  profiles, 
<uv>  normalized  by  u^.  Figure  3.6  shows  that  the  Reynolds  stresses  from  DNS5 
compare  very  well  with  DNS  values  from  Kim  et  al.  (1987).  The  resolved  <uv> 
profile  from  LEST  is  slightly  lower,  as  expected,  but  when  the  SGS  stress  is  added 
{<uv>  +  <ri2>),  the  comparison  is  better. 

Turbulence  intensities 

Wall-normal  distributions  of  turbulence  intensities  show  the  energy  distribution 
throughout  the  channel  for  each  of  the  turbulent  velocity  components.  They  are 
defined  as  the  square  root  of  the  planar  mean  of  the  squares  of  the  fluctuating 
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Figure  3.6;  Reynolds  stress  profiles.  -  DNSS;  .  LES7  (<  uv  >); 

.  LES7  {<uv>  +  <ti2>);  o  Kim  et  al.  (1987). 

velocity  components.  Figure  3.7  shows  the  turbulence  intensities  for  LES7  and 
DNSS  normalized  by  Ur-  The  DNSS  data  compares  very  well  with  the  DNS  data 
of  Kim  et  al.  (1987).  The  LES7  data  has  slightly  higher  values  of  <uu>,  but  lower 
values  of  <vv>  and  <ww>. 

Higher  order  moments 

The  third  moment  of  the  probability  distribution,  the  skewness  factor,  denoted 
by  S,  is  a  measure  of  the  anisotropy  of  the  high  amplitude  velocity  fluctuations. 
For  instance,  a  positive  skewness  factor  indicates  a  preference  for  high  amplitude 
positive  fluctuations  (relative  to  the  local  r.m.s.).  In  this  case,  the  skewness  factors 
have  been  computed  with  the  formula, 

=  ^  (3.1.4) 
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Figure  3.9;  Velocity  kurtosis.  -  DNSS;  .  LES7; -  Kim  et 

al.  (1987).  (a)  u  (b)  v  (c)  w. 

Figure  3.8  shows  very  good  agreement  between  DNSS,  LES7  and  the  DNS  of  Kim  et 
al.  (1987).  The  skewness  factor  for  w  is  very  close  to  zero,  as  required  by  spanwise 
symmetry,  and  is  not  shown. 

The  fourth  moment  of  the  probability  distribution  is  the  kurtosis  or  flatness 
factor,  F,  which  was  computed  directly  from  the  velocity  fields, 

=  (3.1.5) 

The  kurtosis  is  a  measure  of  the  probability  of  high  amplitude  fluctuations.  Fig¬ 
ure  3.9  shows  the  the  kurtosis  distributions  for  the  three  velocity  components.  The 


43 


kurtosis  value  for  a  normal  distribution  [F  =  3),  is  shown  for  comparison.  The 
kurtosis  for  u  is  relatively  near  to  3,  whereas  the  distributions  for  v  and  w  are  very 
large  near  the  wall,  but  converge  to  values  near  3  further  from  the  wall.  The  very 
high  values  near  the  wall  indicate  processes  which  have  very  large  and  infrequent 
excursions  from  the  local  r.m.s..  The  comparison  between  DNSS,  LES7  and  Kim  et 
al.  (1987)  is  very  good. 

Spectra 

Since  one  of  the  objectives  of  this  work  is  to  relate  the  velocity  fields  to  the  wall 
pressures,  in  wavenumber  space,  validation  of  the  velocity  spectra  is  important. 
The  spectra  of  turbulence  has  three  ranges  [see  e.g.,  Bradshaw  (1971)].  First, 
the  “energy-containing”  range  which  contains  the  largest  eddies,  and  is  dependent 
upon  the  geometry,  not  upon  viscosity.  Second,  a  “dissipating”  range  which  con¬ 
tains  the  smallest  eddies  and  is  dominated  by  viscosity;  if  there  is  a  large  enough 
wavenumber  separation  between  the  dissipating  and  energy-containing  ranges  (i.e., 
the  Reynolds  number  is  high  enough),  the  turbulence  in  the  dissipating  range  may 
be  isotropic  and  universal,  be.,  independent  of  the  anisoptropic  turbulence  in  the 
energy-containing  range.  Third,  there  may  exist  an  “inertial  subrange”  between 
the  energy-containing  and  dissipating  ranges;  the  spectra  in  the  inertial  subrange 
decays  as  Velocity  spectra  were  computed  for  many  wall-normal  locations 

for  both  LES7  and  DNSS.  The  spectra  were  averaged  over  approximately  50  real¬ 
izations  and  normalized  by  u^S.  Only  spectra  from  =  5  and  =  150  will  be 
shown  as  these  are  the  same  locations  shown  in  Kim  et  al.  (1987). 

Figures  3.10  through  3.13  show  that  the  streamwise  and  spanwise  spectra 
of  the  turbulent  velocity  fluctuations  compare  very  well  with  those  of  Kim  et 
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al.  (1987).  LES7  compares  very  well  with  DNSS,  although  for  =  150  there  is 
some  attenuation  in  LES7  at  the  highest  spanwise  wavenumbers.  Figure  3.10  shows 
that  for  the  simulations,  the  region  with  a  decay  rate  is  negligible,  which 

is  expected  for  low  Reynolds  number  flows.  At  high  wavenumbers  the  decay- 
rates  increase  with  wavenumber,  indicative  of  well  resolved  simulations  (under¬ 
resolution  would  result  in  a  an  energy  build-up  in  the  highest  wavenumbers).  The 
dynamic  SGS  model  depends  on  resolved  velocity  field  information  between  the 
grid  filter  A  and  the  test  filter  A.  In  the  present  LES,  A  =  2 A,  which  gives 
a  streamwise  wavenumber  cutoff  due  to  A  of  =  6.  Figure  3.10  shows  that 
the  range  6  <  k^S  <  12  is  in  the  dissipating  range  (i.e.  at  wavenumbers  higher 
than  where  k^^/^  decay  rate  is  located).  The  spanwise  spectra  of  u  has  a  peak 
which  corresponds  to  a  wavelength  of  about  100  wall  units  which  is  about  the 
expected  spacing  for  low-  and  high-  speed  streaks.  The  peak  in  the  spanwise 
spectra  of  v  corresponds  to  a  wavelength  of  about  60  wall  units  which  could  be 
due  to  streamwise  vortices.  The  velocity  spectra,  in  terms  of  the  wall  pressure 
source  terms,  will  be  explored  in  depth  in  Chapter  5. 
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Figure  3.10:  Velocity  spectra  at  =  5.  -  DNSS; 

- Kim  et  al.  (1987).  (a)  Streamwise,  (b)  spanwise. 


Figure  3.11;  Spectra  of  udXy^  =  150.  -  DNSS; 

al.  (1987).  (a)  Streamwise,  (b)  spanwise. 


LES7;  -  - 


Kim  et 


Figure  3.12;  Spectra  of  u  at  =  150.  -  DNS5; .  LES7; - Kim  et 

al.  (1987).  (a)  Streamwise,  (b)  spanwise. 
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Figure  3.13;  Spectra  of  w  at  ?/+  =  150.  -  DNS5; .  LES7; - Kim  et 

al.  (1987).  (a)  Streamwise,  (b)  spanwise. 

3.2  Wall-pressure  fields 


In  this  section  we  will  show  the  statistics  and  wavenumber-frequency  spectra  from 
multiple  turbulence  simulations.  We  will  also  describe  how  the  spectra  were  com¬ 
puted. 
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Simulation 

m.s. 

Skew 

Kurtosis 

LES7 

2.05 

-0.06 

5.37 

DNS5 

2.49 

-0.10 

5.23 

Table  II;  Statistics  for  wall  pressure  from  turbulence  simulations.  The  m.s.  values 
have  bee  normalized  by  r^, 

Statistics  of  wall  pressure 

Table  II  is  a  tabulation  of  the  wall  pressure  statistics.  All  the  mean-square  (m.s.) 
values  are  between  2  and  3,  close  to  those  computed  by  Kim  et  al.  (1987)  and 
Kim  (1989).  Handler  (1984)  obtained  a  value  of  11  on  a  coarse  grid  {Ax^  = 
31)  DNS,  while  Moin  and  Kim  (1982)  obtained  4.2  on  a  short  box  =  2ti5) 
LES.  The  values  presented  here  are  much  lower  than  those  obtained  from  exper¬ 
iments  conducted  at  much  higher  Reynolds  numbers.  For  instance,  Farabee  and 
Casarella  (1991)  for  a  boundary  layer  flow  at  =  1535  obtained  m.s.  values  of 
about  9.  Singer  (1996)  obtains  a  value  of  7.5  for  a  simulated  boundary  layer  flow 
with  Rcr  ~  1000.  However,  they  are  consistent  with  the  trend,  shown  by  Farabee 
and  Casarella  (1991),  of  decreasing  m.s.  as  Rsr  decreases. 

The  skewness  is  an  indicator  of  the  asymmetry  of  high  amplitude  fluctuations. 
While  it  has  been  shown  that  the  skewness  of  the  turbulent  pressure  fluctuations 
is  negative  within  the  flow  [see  e.g.^  Kim  et  al.  (1987)  and  Kim  (1989)]  possibly 
due  to  negative  pressure  within  the  vortex  cores,  the  skewness  approaches  very 
small  negative  values  at  the  wall.  It  is  reasonable  that  the  skewness  become  less 
negative  at  the  wall,  because  the  probability  of  high  amplitude  positive  pressure 
fluctuations,  corresponding  to  sweep  motions,  increases  as  the  wall  is  approached. 
Schewe  (1983)  showed  that  as  transducer  size  increases  (say  for  >  100),  skew¬ 
ness  approaches  the  Gaussian  limit  of  0;  as  transducer  size  decreases  skewness 
approaches  a  value  of  approximately  —0.2.  As  shown  in  Table  II  the  skewness  for 
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the  simulations  are  small  negative  numbers. 

The  kurtosis  shows  that  the  wall  pressure  is  a  highly  intermittent  function 
where  a  large  proportion  (large  compared  with  a  Gaussian  distribution)  of  the 
r.m.s.  pressure  comes  from  high  amplitude  fluctuations.  Schewe  (1983)  showed 
values  of  kurtosis  approaching  5  as  the  pressure  transducer  size  decreased,  while 
approaching  the  Gaussian  limit  of  3  as  the  size  increased.  However,  one  may  also 
expect  kurtosis  to  increase  with  increasing  grid  size.  Piomelli  (1993)  explains  that 
high  kurtosis  for  velocities  may  be  the  result  of  not  resolving  the  smaller  scales 
of  turbulence.  These  scales  tend  to  be  more  isotropic  than  the  larger  scales,  and 
if  resolved  would  decrease  the  kurtosis.  The  kurtosis  values  of  wall  pressure  are 
around  5  for  all  the  simulations,  which,  based  on  the  values  of  Ax"'',  are  higher 
than  Schewe’s  trend.  However,  the  kurtosis  from  the  simulations  do  not  show  any 
dependency  on  streamwise  grid  size. 

Spectra 

Computation  of  spectra 

The  wavevector-frequency  spectra  shown  in  this  Section  were  computed  by  the 
following  methodology;  the  data  were  divided  into  Nt  =  384  point  segments  span¬ 
ning  the  period,  T,  with  50  percent  overlap.  This  gave  a  reasonable  compromise 
between  number  of  spectra  realizations  and  frequency  resolution.  A  Hanning  win¬ 
dow  was  used  on  the  time  segments,  while  no  windowing  was  necessary  in  the  x 
and  z  directions  in  which  the  flow  is  periodic.  The  Hanning  window  was  applied  in 
time  to  the  total,  resolved  pressure,  p,  normalized  by  the  long  time  average  value 
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(3.2.6) 


of  T„ 


1 


2tiV 


p  (xi,Zj,  ti)  =  p  {Xi,  Zj,ti)  2  I  ^  ^ 


where  i  =  1 ...  N^,  j  =  1 ...  N^,  and  I  =  1 ...  N^.  A  three-dimensional  Fourier 
wavenumber-frequency  transform  of  each  time  segment  was  taken  and  from  that, 
an  averaged  spectral  density  was  obtained  from  the  mean  of  the  spectral  densities 
from  the  m  overlapping  segments 


1  m 

(iAk^jAk^JAuj)  =  —Y.(PP*)s-  (3.2.7) 

Here  p  denote  the  Fourier  coefficients  of  the  pressure  and  p*  are  their  complex 
conjugates.  The  spectral  density  was  computed  for  both  top  and  bottom  of  the 
channel,  resulting  in  2m  time  segments 


The  three-dimensional  spectral  density  should  satisfy  Parseval’s  identity 

T  Lz  Lx  OO  OO  OO 

ITITt  I  f  J p^  =  j  j  j  dio  dk^  dk^  (3.2.8) 

^  ^  0  0  0  — OO  — OO  — OO 

which,  in  discretized  form,  is  given  by 

1  Nt  Nz  Nt  Nz  Nrc 

.r  .r  .r  P^  t/ )  =  AlO  Ak,  ^  ^  ^  {iAkx,  jAkz,  lAuj) 

(3.2.9) 

where  the  left  hand  side  is  the  mean-square  (m.s.)  pressure.  Since  the  pressure  has 
been  windowed  in  time,  its  spectral  density  does  not  satisfy  (3.2.9)  exactly;  d*  is 
then  normalized  to  satisfy  (3.2.9).  The  two-dimensional  kx  —  ui  spectra  is  given  by 
the  integral  over  all  the  spanwise  wavenumber  modes,  which,  in  discretized  form, 
results  in  the  summation  over  all  the  modes  multiplied  by 

^{iAkx.,lAix))  =  AkzYj^-  (3.2.10) 

i=i 
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The  two-dimensional  spectral  density  extends  over  I  =  1 . . .  Nt  frequency  bins  and 
i  =  1 . . .  Nx  wavenumber  bins  and  has  the  symmetry  property 

lAiv)  =  $  [{Nx  —  i  +  2)Akx,  {Nt  —  I  +  2)Alj]  (3.2.11) 

In  the  kx  —  (jJ  plots  the  wavenumber  bins  Nx/2  +  2  <  i  <  Nx  have  been  translated 
to  the  bins  —Nx/2  -|-  1  <  i  <  —  1;  with  this  convention,  a  wave  that  is  convecting 
downstream  will  have  either  u)  or  kx  negative.  In  the  two-dimensional  spectra 
plots,  the  convective  ridge,  representing  energy  convecting  downstream,  is  shown 
in  the  negative  kx,  positive  u)  quadrant. 

DNSl  is  the  database  used  Choi  and  Moin  (1990),  which  was  reduced  by  Bruce 
Abraham  of  Naval  Underwater  Systems  Command  and  presented  here  as  well  as 
in  Chang  et  al.  (1994).  The  DNSl  spectra  were  computed  using  the  MATLAB* 
software  package.  The  data  were  divided  into  four  512  point  time  segments  with 
no  overlap  and  the  two-dimensional  spectral  density  for  each  segment  computed. 
This  two-dimensional  kx  —  to  spectral  density  was  then  averaged  over  the  m  =  4 
segments  and  over  every  other  spanwise  location: 

O  Nz  r  -I  m 

^{iAkx,lAuj)  =  —  ^  —^^{iAkx,jAk^,lAu},r)  .  (3.2.12) 

j  even 

Two-dimensional  wall  pressure  spectra 

One  of  the  reasons  for  running  the  simulations  was  to  study  the  fidelity  of  the 
wavenumber-frequency  characteristics  of  the  wall  pressure  fields,  particularly  in 
the  subconvective  wavenumber  rangeb  Table  III  lists  the  number  of  grid  points, 
*  MATLAB  is  a  registered  trademark  of  The  MathWorks,  Inc. 

^The  subconvective  wavenumber  range  is  the  region  of  the  wavenumber-frequency  plane  where 
kx  kxc  where  kxc  =  i^/Uc,  is  the  location  of  the  convective  ridge;  Uc  is  the  convection  velocity 
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Simulation 

Grid  Points 
N^x  NyX 

Box  Size 
%  X  %  X  % 

^  d  d 

Dealiasing 

LESl 

48  X  64  X  64 

none 

LES2 

48  X  64  X  64 

X 

LES3 

48  X  64  X  64 

X,  z 

LES7 

72  X  64  X  64 

X,  z 

DNSl 

128x129x128 

X,  z 

DNS2 

64  X  64  X  64 

none 

DNS3 

192  X  96  X  96 

none 

DNS5 

192  X  96  X  96 

X,  z 

Table  III:  Simulation  designations  and  details 
the  box  sizes  and  dealiasing  used  in  the  simulations.  To  understand  why  the 
various  simulations  were  run,  the  two-dimensional  wavenumber-frequency  spectra 
of  wall  pressure  will  be  shown,  since  it  was  found  that  numerical  issues  (namely, 
aliasing  errors  in  the  pseudo-spectral  method  as  well  as  grid  resolution)  affected 
the  behavior  of  the  spectra,  particularly  in  the  subconvective  wavenumber  range. 


LESl  and  LES2  were  run  on  a  streamwise  box  size,  whereas  LES7  was 
run  on  a  streamwise  box  size.  Eigure  3.14  shows  plots  of  the  wavenumber- 
frequency  spectra  of  wall  pressure.  The  wavenumber-frequency  spectra  is  displayed 
as  contours  of 


10  log 


<l>(zAfci,  1/S.u) 


(3.2.13) 


Tl6‘^/Ur 

Eigure  3.14(a)  the  wavenumber-frequency  spectra  from  LESl  has  energy  in  the 
positive-A:^;  quadrant  which  was  shown  by  Chang  et  al.  (1994)  to  be  due  to  aliased 
energy.  LES2  and  LES7  were  run  with  dealiasing  in  x,  and  x  —  z,  respectively. 
As  can  be  seen  from  Eigure  3.14(b)  and  Figure  3.15,  the  energy  in  the  positive-A:^ 
quadrant  is  significantly  reduced  by  dealiasing. 


The  wavenumber-frequency  spectra  for  DNSl  are  shown  in  Figure  3.16.  The 


and  uj  denotes  the  radian  frequency. 
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Figure  3.14:  kx  —  oo  contour  plots  of  wall  pressure  spectra  from  (a)  LESl  and  (b) 
LES2. 

high  energy  build-up  around  kxS  =  0  and  the  high  subconvective-wavenumber 
energy  levels  could  not  be  explained  by  Choi  and  Moin  (1990).  DNS2  was  run  on 
a  very  short  box  (271(5)  with  a  high  grid  density  to  see  the  effect  of  grid  resolution 
on  the  wavenumber-frequency  characteristics.  As  can  be  seen  from  Figure  3.17, 
DNS2  also  presents  an  energy  build-up  in  the  positive- /la,  quadrant,  and  high  energy 
levels  in  the  subconvective  wavenumber  range.  DNS3  was  run  on  a  QnS  box  without 
dealiasing;  as  Figure  3.18  shows,  it  has  energy  in  the  positive-/Oa,  quadrant,  due  to 
aliasing,  as  well  as  an  energy  buildup  in  the  subconvective  wavenumber  range. 
DNSS  was  dealiased  in  both  x  and  z,  and  as  shown  in  Figure  3.19,  it  has  no  energy 
in  the  positive-Zi^  quadrant  (except  for  that  related  to  the  convective  ridge)  and 
shows  no  energy  buildup  in  the  subconvective  wavenumber  range. 
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Figure  3.18;  kx  —  lo  contour  plot  of  wall  pressure  spectra  from  DNS3 


Chapter  4 


Partial  Wall-Pressures 


The  objective  of  this  chapter  is  to  determine  which  velocity  source  terms  are  most 
important  for  the  generation  of  wall  pressure,  and  where  in  the  channel  these 
sources  are  acting.  In  order  to  do  this,  partial  wall  pressures,  i.e.,  the  contributions 
due  to  individual  source  terms  at  various  horizontal  layers  of  the  channel,  have  been 
computed  and  are  compared  in  planar  wavenumber  space.  Partial  pressures  are 
the  product  of  the  wall  pressure  Green’s  function  and  the  individual  source  terms, 
integrated  over  vertical  regions,  corresponding  to  various  regions  of  the  channel. 
Comparisons  of  the  partial  pressures  show  which  sources  are  most  significant,  while 
partial  pressures  are  combined  and  compared  to  the  total  pressure  to  determine 
which  subset  of  sources  is  necessary  to  reconstitute  the  total  pressure. 

Section  4.1  describes  how  the  partial  pressures  are  computed;  Section  4.2  shows 
the  relative  importance  of  the  Stokes’  pressures;  Section  4.3  compares  the  total, 
TT  and  MS  partial  pressures;  Section  4.4  compares  the  partial  pressures  due  to 
individual  terms  and  regions;  Section  4.5  summarizes  the  hndings. 


57 


Region 

Limits 

Description 

1 

0  <  ?/+  <  5 

Viscous  shear-layer 

2 

5  <  ?/+  <  30 

Buffer  layer 

3 

30  <  ?/+  <  180 

Logarithmic  region 

4 

180  <y+  <  360 

Upper  channel 

Table  I;  Regions  of  the  channel. 

4.1  Computation  of  partial  pressures 

Since  the  wall  pressure  depends  on  sources  distributed  throughout  the  entire  flow 
domain  it  is  difficult  to  ascertain  where  in  the  boundary  layer  the  dominant  sources 
are  located.  One  advantage  of  turbulence  simulations  is  that  they  generate  three- 
dimensional  velocity  fields  from  which  all  the  source  terms  can  be  computed.  In 
this  investigation  the  product  of  the  sources  by  the  Green’s  function,  which  is  the 
true  measure  of  the  influence  of  the  sources  on  the  wall  pressure,  has  been  com¬ 
puted.  The  channel  was  divided  into  four  regions  (Table  I)  roughly  corresponding 
to  the  viscous  shear-layer,  buffer  and  logarithmic  regions,  and  the  upper  part  of 
the  channel.  The  contribution  of  the  source  terms  was  integrated  over  each  region. 
It  is  assumed  that  the  partial  pressures  are  evaluated  at  the  lower  channel  wall, 
even  though  symmetry  has  been  used  in  the  data  reduction  and  partial  pressures 
have  been  averaged  over  both  the  upper  and  lower  channel  walls.  The  lower  wall 
will  be  referred  to  as  the  “near-wall,”  the  upper  wall,  the  “far-wall.” 

The  source  terms,  Tij{x,y,  z),  were  computed  from  110  realizations  of  DNSS 
data,  in  three-dimensional  physical  space,  using  the  2/3  rule  to  eliminate  aliasing 
errors.  Nine  terms  were  computed;  the  six  TT  terms,  their  sum,  the  MS  term  and 
the  total.  They  were  then  Fourier-transformed  to  yield  their  complex  coefficients 


in  planar  wavenumber  space,  Tjj(K,  y).  The  integral  of  the  product  of  the  source 
terms  and  Green’s  function 

MK,Rr)  =  fijiK,7])G{K,y  =  -l,y)dy.  (4.1.1) 

J-ni{r) 

(where  r^i(r)  and  r]2{r)  are  the  lower  and  upper  integration  limits  for  each  region, 
and  r  denotes  the  vertical  regions  described  in  Table  I)  was  then  computed  using 
Chebychev  integration.  The  complex  partial  pressure  spectrum  due  to  the  source 
term  Rj  from  region  r  is  denoted  by  p(K,  Rr).  The  complex  spectra  for  the  various 
terms  and  regions  could  then  be  combined;  e.g.,  the  partial  pressure  due  to  T12 
from  regions  2,  3  and  4  is 

P^,pK,mU)=pq(K,R2)+p^(K,m)+pS(Vi,RA).  (4.1.2) 

The  total  TT  partial  pressure  is  the  summation  of  all  the  TT  partial  pressures, 

r^(K,i?r)  =  ^^p5^(K,7?r),  (4.1.3) 

and  the  total  partial  pressure  is 

^“*(K,/?r)  =/^^(K,i?r)+p^^(K,/?r).  (4.1.4) 

The  total  pressure  is  |?°*(K,  /?1234).  The  magnitudes  of  the  complex  spectra  were 
computed  by 

n(K, /?r)  =  |1*^  (K,  i?r)|ljj(K,  i?r)  (no  summation)  (4.1.5) 

(where  f*  is  the  complex  conjugate  of  /),  and  averaged  over  multiple  realizations. 
It  should  be  noted  that  the  partial  pressures  are  defined  by  the  integral  over  all 
frequencies, 

/OO 

n{K, uj,Rr)duj  (4.1.6) 

-OO 
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where  n(K  ,a;,  Rr)  is  the  partial  pressure  in  four-space.  The  one-dimensional  spec¬ 
tra  were  obtained  by  the  integrations 

/OO 

n(K,  (4.1.7) 

-OO 

and 

/OO 

U{¥.,Rr)dK  (4.1.8) 

-OO 

The  one- dimensional  spectra  were  normalized  by  t‘^5. 


4.2  Stokes  Pressure 


The  partial  pressure  due  to  the  boundary  terms  (2.6.17)  is  called  the  Stokes  pres¬ 
sure  [Kim  (1989)].  Kraichnan  (1956)  estimated  that  the  Stokes  pressure  should  be 
much  smaller  than  the  total  (due  to  velocity  field  sources)  pressure.  Kim  (1989), 
using  DNS  data  from  a  channel  flow,  computed  the  r.m.s.  of  the  Stokes  wall  pres¬ 
sure  and  found  it  to  be  about  1/10  of  the  total  pressure.  In  his  analytical  model, 
Chase  (1980)  assumed  it  was  small  and  neglected  it.  In  this  subsection  we  will 
show  the  wavenumber  distribution  of  the  Stokes  pressure  and  compare  it  to  the 
total  pressure. 


The  pressure  gradients  in  (2.6.17)  can  be  expressed  in  terms  of  velocities  as  a 
result  of  the  wall  normal  momentum  equation  evaluated  at  the  walls. 


dp 

dy 


y=±l 


1  dd'v 


Re  dy'^ 


y=l 


(4.2.9) 


In  planar  wavespace,  the  far-wall  boundary  term,  evaluated  at  the  near-wall,  be¬ 
comes 


B^iK,y  =  -l)  = 


1 


1  d^v 


K  sinh(2Kr)  Re  drf 


(K,77  =  1) 


(4.2.10) 
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Figure  4.1;  Partial  near- wall  wall  pressure  due  to  boundary  condition  terms  and  to¬ 
tal  pressure  with  and  without  boundary  condition  terms;  :  near- wall  boundary 

condition  term;  A  :  far-wall  boundary  condition  term; - :  ptot  with  boundary 

condition  terms; .  :  ptot  without  boundary  condition  terms;  (a)  kx  (b)  kz- 

and  the  near-wall  boundary  term,  evaluated  at  the  near-wall,  becomes 


M  Ihi  k5 


B^{K,y  =  -l) 


1  1 

K  tanh(2A^)  Re  ^ 


1). 


(4.2.11) 


The  magnitudes  of  the  complex  spectra  (4.2.10)  and  (4.2.11),  were  taken  and 
then  integrated  appropriately  to  obtain  either  the  streamwise  or  spanwise  one¬ 
dimensional  spectra.  Figure  4.1  shows  the  one-dimensional  spectra  for  the  near- 
and  far-wall  boundary  terms  computed  from  11  widely  spaced  (in  time)  realizations 
of  DNSS.  The  streamwise  spectra  show  the  near-wall  term  is  about  15  —  20  dB 
less  than  the  total  pressure.  The  difference  between  the  total  pressure  with  and 
without  the  boundary  terms  is  less  than  1  dB  throughout  the  entire  streamwise 
wavenumber  range.  The  far-wall  boundary  term  is  about  18  dB  lower  than  the 
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total  pressure,  at  the  lowest  wavenumber,  and  then  falls  off  rapidly.  The  spanwise 
spectra  shows  that  the  near-wall  term  comes  to  within  2^3  dB  of  the  total 
pressure  at  the  high  wavenumbers.  The  difference  between  the  total  pressure  with 
and  without  the  boundary  term  is  as  great  as  2  dB.  The  conclusion  that  can 
be  drawn  from  this  is  that,  for  the  wall  pressure  at  the  near-wall,  the  far-wall 
boundary  term  is  insignificant,  while  the  near-wall  boundary  term  has  a  small,  but 
discernible  affect  on  the  near-wall  wall  pressure  spectra.  For  all  practical  purposes, 
neglecting  the  Stokes  pressure  is  a  reasonable  assumption. 

4.3  Total,  MS  and  TT  partial  pressures 

All  regions 

It  was  first  estimated  by  Kraichnan  (1956)  that  the  mean-square  of  the  MS  pressure 
should  be  approximately  68  times  larger  than  the  TT  pressure.  Using  DNS  data 
however,  Kim  (1989)  found  that  the  one-dimensional  spectra  of  the  MS  and  TT 
pressures  were  the  same  order  of  magnitude  throughout  the  wavenumber  range. 
One  reason  that  Kraichnan’s  (1956)  estimate  does  not  compare  well  with  more 
recent  data  may  be  because  it  was  dependent  upon  estimates  of  fluctuating  velocity 
correlation  lengths  obtained  for  homogenous,  isotropic  turbulence. 

Figure  4.2  shows  the  total,  MS  and  TT  pressures  summed  over  all  regions, 
'Kij{kx^  Rl'2,‘iA)  and  7rjj(A;2,  i?1234).  Consistent  with  the  results  of  Kim  (1989),  the 
MS  and  TT  terms  are  of  the  same  order  of  magnitude  throughout  the  wavenum¬ 
ber  range.  The  MS  term  is  larger  than  the  TT  term  throughout  the  streamwise 
wavenumber  range  while  in  the  spanwise  spectra,  for  k^S  >  25,  the  TT  term  be- 
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comes  larger. 

The  low- wavenumber  range  of  the  streamwise  spectra  is  fairly  flat,  with  a  peak 
at  kx5  ~  3.  Farabee  and  Casarella  (1991)  used  u-r,  and  5  (where  5  is  the 
boundary  layer  thickness)  to  normalize  the  frequency  spectra,  and  found  that, 
for  a  range  of  Reynolds  numbers,  the  spectra  collapsed  in  a  mid-frequency  range. 
The  peak  of  their  spectra  was  at  ujUt/S  ^  50.  Assuming  a  convection  velocity 
Uc  ~  l?)Ur  (based  on  the  location  of  the  convective  ridge  in  the  —  to  plane)  the 
spectral  peak  they  obtained  corresponds  to  kxS  ~  4,  in  good  agreement  with  the 
data  presented  here. 

The  total,  MS  and  TT  pressures  do  not  appear  to  have  significant  regions  of 
k^^  or  k^^  behavior  predicted  by  Bradshaw  (1967)  for  the  overlap  region.  This  is 
due  to  Reynolds  number  effects;  Farabee  and  Casarella  (1991)  define  the  overlap 
region  as 

100  <  to5/ur  <  O.mRer  (4.3.12) 

which  for  Rcr  =  180  and  Uc  ~  ISu^,  gives  a  lower  limit  kxS  ~  8  and  an  upper 
limit  kxS  ~  4. 

Both  the  streamwise  and  spanwise  spectra  have  wavenumber  ranges  with  a  k^^ 
decay  rate.  This  range  has  been  shown  to  scale  with  inner  variables  [Farabee  and 
Casarella  (1991),  Keith  et  al.  (1992),  Olivero-Bally  (1993),  Chang  et  al.  (1994)] 
indicating  that  it  is  due  to  sources  in  a  region  that  scales  accordingly,  such  as 
the  viscous  shear-layer.  The  MS  spanwise  spectra  has  a  significant  range  of  k^^ 
whereas  the  TT  spanwise  spectra  only  passes  through  a  rolloff.  This  will  be 
investigated  further  in  the  next  subsection. 
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Figure  4.2:  Total,  MS  and  TT  one-dimensional  pressure  spectra.  — 
. .  (a)  Streamwise,  (b)  spanwise. 

Contribution  of  the  regions 

To  determine  where  in  the  boundary  layer  the  MS  and  TT  pressures  are  generated, 
we  will  now  examine  the  partial  pressures  from  the  four  regions. 

Total  partial  pressures 

The  spectra  of  the  contribution  to  the  total  pressure  from  each  of  the  four  regions, 
7r*°*(i?r),  and  for  all  four  regions  combined,  7r‘"‘(i?,1234),  are  shown  in  Figure  4.3. 
For  kxS  >  5  and  kz5  <  50  the  buffer  layer  accounts  for  most  of  the  total  pressure 
spectrum;  for  k^S  >  55,  however,  the  viscous  shear-layer  becomes  the  most  im¬ 
portant  contribution.  The  high-wavenumber  behavior  of  the  viscous  shear-layer  is 
consistent  with  the  notion  that  the  near-wall  regions  are  responsible  for  the  higher- 
wavenumber  pressure,  since  the  high- wavenumber  contributions  to  the  pressure 
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Figure  4.3:  One-dimensional  spectra  of  the  total  pressure  for  the  various  re¬ 
gions.  All  regions  (R1234); .  viscous  shear-layer  (Rl); .  buffer 

layer (R2); - logarithmic  region(R3); - upper  channel(R4).  (a)  Stream- 

wise,  (b)  spanwise. 

from  regions  further  from  the  wall  are  strongly  attenuated  by  the  Green’s  func¬ 
tion.  At  low  wavenumbers  <  5,  <  5)  the  partial  pressures  show  very  little 

resemblance  to  7r*°*(/?1234).  There,  the  total  source  term  shown  in  Figure  4.4 
at  several  y  locations,  is  nearly  flat,  and  (i?r)  must  have  the  same  slope  as  the 
Green’s  function,  as  shown  in  Figure  4.3. 

MS  partial  pressures 

In  Figure  4.5  the  spectra  of  the  MS  pressure  are  shown  for  each  of  the  four  regions. 
Except  at  the  lowest  wavenumbers,  TT^^^^kx)  is  due  primarily  to  buffer  layer  sources. 
This  is  also  true  for  7r^^^{kz),  except  at  the  highest  wavenumbers  {k^S  >  50),  where 
the  viscous  shear-layer  contributions  become  significant.  It  has  been  conjectured 
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Figure  4.4:  Streamwise  wavenumber  spectra  of  total  source  terms  at  selected  loca¬ 
tions.  -  =  5; .  y'^  =  11; .  y^  =  21;  A  =  49;  +  y'^  =  145. 


Figure  4.5;  One-dimensional  spectra  of  the  MS  pressure  for  the  various  regions. 

All  regions  (R1234); .  viscous  shear-layer  (Rl);  .  buffer  layer 

(R2); - logarithmic  region  (R3); - upper  channel  (R4).  (a)  Streamwise, 

(b)  spanwise. 
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{e.g.,  by  Blake  (1986),  p.  518)  that  at  the  low  wavenumbers  the  MS  sources 
follow  a  k'^  behavior.  Figure  4.5(a)  shows  that  in  the  low  wavenumbers  the  partial 
pressure  from  the  buffer  layer  has  a  k'^  slope.  The  figure  also  shows  that  the 
buffer  layer  also  has  a  very  short  region  with  a  k~^  slope,  which  is  consistent  with 
Bradshaw’s  (1967)  prediction.  The  partial  pressure  from  the  viscous  shear-layer 
has  a  range  of  k~^  behavior,  which  conhrms  that  the  approximations  involved  in 
MS  source  term  modelling  in  the  viscous  shear-layer  [Blake  (1986)]  are  essentially 
correct.  However,  as  can  be  seen,  in  the  k~^  wavenumber  range,  the  sources  from 
the  buffer  layer  are  dominant,  rather  than  those  from  the  viscous  shear-layer.  It 
can  also  be  seen,  in  Figure  4.5(b),  that  the  spanwise  partial  pressure  spectra  due 
to  the  buffer  layer  has  wavenumber  ranges  that  decay  as  k^^  and  k^^. 

Now  that  we  know  which  regions  dominate  the  MS  partial  pressure  it  is  of 
interest  to  see  which  combinations  of  regions  are  necessary  to  obtain  the  total  MS 
partial  pressure.  This  is  particularly  of  interest  for  modelling  purposes,  since  it 
would  give  indications  as  to  which  minimum  subset  of  regions  is  necessary  to  model 
the  MS  pressure.  Figure  4.6  shows  the  MS  partial  pressures  for  combinations  of  the 
regions.  For  discussion  purposes,  wavenumber  ranges  are  defined  in  Table  II.  In  the 
intermediate  and  high  streamwise  and  spanwise  wavenumbers,  contributions  from 
the  viscous  shear-layer  and  buffer  layer  are  necessary  to  reconstitute  the  MS  partial 
pressure.  This  confirms  the  dominance  of  the  buffer  layer  and  refutes  the  notion 
that  only  the  viscous  shear-layer  sources  are  important  at  the  high  wavenumbers. 
In  the  lowest  and  low  wavenumber  range  only  the  buffer  layer  and  logarithmic 
region  sources  are  necessary. 

TT  partial  pressures 
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Figure  4.6;  One-dimensional  spectra  of  the  MS  pressure  for  combinations  of  re¬ 
gions.  All  regions  (R1234);  A  viscous  shear-layer,  buffer  layer  and  logarith¬ 
mic  region(R123);  -|-  viscous  shear-layer  and  buffer  layer(R12)  o  buffer  layer  and 
logarithmic  region(R23).  (a)  Streamwise,  (b)  spanwise. 


Spectra 

Wavenumber  range 

Lowest 

Low 

Intermediate 

High 

kxS  <  1 

1  <  kxS,  kz5  <  5 

5  <  kx5^  kzS  <  30 

30  <  kz6  <  70 

^MS 

2+3 

2+3 

1+2 

1+2 

^TT 

1+2+3+4 

1+2+3 

1+2 

1+2 

Table  II;  Regions  of  channel  which  dominate  the  MS  and  TT  spectra.  1;  viscous 
shear-layer;  2;  buffer  layer;  3;  logarithmic  region;  4;  upper  channel. 
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Figure  4.7:  One-dimensional  spectra  of  the  TT  pressure  for  the  various  regions. 

All  regions  (R1234);  .  viscous  shear-layer (Rl);  .  buffer  layer 

(R2); - logarithmic  region  (R3); - upper  channel  (R4).  (a)  Streamwise, 

(b)  spanwise. 

Figure  4.7  shows  the  spectra  of  the  partial  contributions  to  the  TT  pressure  for 
each  of  the  four  regions.  In  this  case  the  individual  curves  have  a  significantly 
different  behavior  than  their  sum.  The  buffer  layer  gives  the  largest  contribution, 
except  at  the  highest  spanwise  wavenumbers,  where  the  viscous  shear-layer  be¬ 
comes  important.  Since,  at  low  wavenumbers,  the  sum  of  the  contributions  to  the 
total  pressure  from  the  various  regions  results  in  spectral  values  much  lower  than 
those  of  each  region,  the  contributions  of  the  individual  regions  must  cancel  each 
other  out,  indicating  that  long- wavelength  turbulent  structures  may  have  coher¬ 
ence  over  the  entire  channel.  The  high  levels  of  the  total  pressure  spectra  at  low 
wavenumbers  shown  previously  in  Figure  4.3  are  due  to  the  TT  contribution,  since 
the  MS  pressure  spectra  were  previously  shown  to  be  flat  at  low  wavenumbers. 
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Figure  4.8:  Wall-normal  distribution  of  77)  .  -  =  1/3,  kz5  =  0; 

.  kj  =  1/3,  M  =  3/2;  —  kj  =  1,  M  =  9/2; - =  3,  k,5  =  12. 

The  high  levels  of  the  total-  and  TT-pressure  spectra  from  the  upper  channel 
that  are  seen  in  Figures  4.3(a)  and  4.7(a)  can  be  explained  by  examining  the 
integrand  of  (4.1.1),  =  h*^hij,  where  hij  is  the  product  of  the 

source  term  and  the  Green’s  function, 

(4.3.13) 

Figure  4.8  shows  for  various  values  of  K.  The  lowest  K  combination, 

which  represents  a  planar  wave  with  wavelength  equal  to  the  streamwise  box  size, 
presents  very  little  attenuation  across  the  channel,  indicating  that  the  pressure  at 
the  lowest  wavenumbers  at  one  wall  is  affected  by  sources  from  the  buffer  layer  on 
the  opposite  one,  he.,  that  the  Green’s  function  is  a  low-pass  filter  that  is  almost 
transparent  to  the  low-wavenumber  sources.  This  does  not  explain,  however,  why 
the  TT  partial  pressures  from  all  the  regions  cancel  each  other  out. 

Figure  4.9  shows  the  TT  partial  pressure  spectra  obtained  by  combinations 
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Figure  4.9:  One-dimensional  spectra  of  the  TT  pressure  due  to  combinations  of 
regions.  All  regions  (R1234);  A  viscous  shear-layer,  buffer  layer  and  loga¬ 

rithmic  region  (R123);  -f  viscous  shear-layer  and  buffer  layer  (R12)  o  buffer  layer 
and  logarithmic  region  (R23).  (a)  Streamwise,  (b)  spanwise. 

of  the  various  regions.  The  sum  of  the  partial  pressures  due  to  viscous  shear- 

layer  and  buffer  layer,  7r^^(i?12),  matches  7r^^(i?1234)  for  intermediate  and  high 

wavenumbers  {kxS,kzS  >  5).  In  the  low- wavenumber  range,  however,  the  sources 

from  the  logarithmic  region  become  important,  as  shown  by  the  agreement  between 

n'^'^{R123)  and  7r"^"^(/?1234).  The  present  results  indicate  that  the  viscous  shear- 

layer  gives  a  very  signihcant  contribution  to  the  turbulent-pressure  spectra:  the 

partial-pressure  spectrum  due  to  the  buffer  and  logarithmic  layers  alone,  7r^^(i?23), 

diverges  from  7r^^(i?1234)  and  7r^^(i?123)  even  at  low  wavenumbers.  All  regions 

give  significant  contributions  to  the  lowest  wavenumbers,  again  due  to  the  small 

amount  of  attenuation  provided  by  the  Green’s  function.  The  partial-pressure 

spectra  due  to  each  region  are  higher  than  those  due  to  their  combination,  as  also 

observed  before.  Table  II  summarizes  the  regions  that  were  found  to  affect  most 
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the  spectra  of  the  MS  and  TT  pressure.  The  wavenumber  ranges  dehned  in  the 
Table  will  be  used  throughout  the  rest  of  this  dissertation. 


4.4  Contribution  of  individual  TT  partial  pres¬ 
sures 

In  this  section  we  will  investigate  which  of  the  TT  components  affects  most  the  TT 
partial  pressure  spectra,  by  first  examining  the  terms  singly,  then  by  combining 
them  in  various  ways  and  comparing  them  to  the  total-  or  TT-pressure  spectra. 

Figure  4.10  shows  the  spectra  of  the  most  significant  TT  terms  and  regions. 
Similar  to  what  was  observed  before,  the  spectra  of  the  individual  TT  terms  have 
higher  low-wavenumber  content  (in  the  streamwise  direction)  than  their  combina¬ 
tion.  In  the  lowest-  and  low- wavenumber  ranges  7r2^{R2),  7r^^(7?3)  and  7rf^{R2) 
are  the  dominant  terms;  all  the  TT  terms,  particularly  those  due  to  the  logarithmic 
and  buffer  regions,  have  similar  behavior  and  are  higher  than  ti'^'^^RA)  at  the  low¬ 
est  wavenumbers.  In  the  range  10  <  <  20,  (/?2)  and  Trfg^ (i?2)  are  the  most 

significant  while  for  20  <  <  30,  three  terms,  7r^^(i?2),  t{(^{R2),  and  7rJ7’(i?2) 

dominate,  n'23  {R2)  is  the  largest  term  for  the  spanwise  spectra  as  well,  except  at 
the  highest  wavenumbers  {k^S  >  60),  where  contributions  from  the  viscous  shear- 
layer  become  important.  Figure  4.11  shows  that  'k'^{R1)^  ti'^^{R1)  and  7r^^(i?l) 
are  the  dominant  terms  in  the  highest  kz  wavenumber  range.  It  is  interesting  to 
observe  that  the  diagonal  terms  7r^^(i?l)  and  {Rl)  decay  like  whereas  the 
off-diagonal  terms,  {R\)  and  7r^^(i?l),  decay  at  a  steeper  rate,  roughly 

Once  the  dominant  terms  and  regions  were  identified,  they  were  combined  to 
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Figure  4.10;  One-dimensional  partial  pressure  spectra  due  to  dominant  TT  source 

terms  and  regions.  7r^^(i?1234);  o  {Rl)\ - tt'2^{R2)]  x  {R3)] 

A  7rJ'.2  {R2); - 7iJ^{R2).  (a)  Streamwise,  (b)  spanwise. 


kz5 

Figure  4.11;  Partial  pressures  for  sources  in  the  viscous  shear-layer.  .  : 

A  :  7ri7(i?l); .  :  nJi'iRiy,  +  : 

determine  whether  subsets  of  terms  or  regions  could  re-constitute  the  total  pres¬ 
sure.  The  combinations  are  summarized  in  Table  III.  Combination  4,  C4,  which 
combines  {R2)  and  {R2)  is  shown  in  Figure  4.12.  The  addition 

of  n^^iRl)  (C17),  which  is  the  dominant  term  at  higher  wavenumbers  (see  Fig¬ 
ure  4.5)  gives  significant  improvement.  Additional  terms,  such  as  7r^^^{R3)  (C18) 
and  7123  (i?3)  and  7r[^ {R3)  (C19)  did  not  give  further  improvements.  The  spanwise 
spectra  extend  to  much  higher  wavenumbers  than  the  streamwise  ones;  to  obtain 
agreement  for  the  highest  wavenumbers  it  was  necessary  to  add  the  contribution 
of  7r33^(i?123).  (In  Figure  4.11,  (Rl)  was  shown  to  be  the  dominant  partial 
pressure  above  kz5  =  60;  in  the  high  wavenumber  range  {R23)  is  much  smaller 
than  713^ (Rl)  and  the  fact  that  7r|^^(i?123)  was  added  instead  of  7rJ^{Rl)  should 
not  make  an  appreciable  difference  to  our  conclusions). 

While  the  dominant  contributions  to  the  high- wavenumber  regions  of  the  pres¬ 
sure  spectra  could  be  identified  and  isolated,  no  such  simplification  could  be  made 
for  the  low- wavenumber  range.  Numerous  combinations  were  tried;  Figure  4.13 
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Combination 

Terms 

Cl 

pU(i?2)+p?J(i?2) 

C2 

pW{m)+pWm 

C3 

pU(i?23)+p?J(i?23) 

C4 

p^^^R2)+p^i{R2)+pJ^{R2) 

C17 

p^^^R12)+p^i'{R2)  +pfi'iR2) 

CIS 

p^^^Rl23)+p'^i'iR2)+pJi'iR2) 

C19 

p^^^(i?123)  +pU(i?23)  +pWiR‘2S) 

C26 

p^^^{Rl23)+p^i'iRl23)+p{i'{Rl23) 

C27 

C17  +  pW{R123) 

C20 

p^'^{R123)  -pJ[{R123) 

C21 

p^^{R123)-p^nRm 

C22 

p^^{R123)-pJl{R123) 

C23 

p^^{R123)  —  pf^(i?123) 

C24 

p^^{R123)  —  pj^  {R123) 

C25 

p'^^{R123)-p^^{R123) 

Table  III;  Combinations  of  terms  used  in  attempts  to  reconstitute  the  total  wall 
pressure  spectra. 

shows  that  even  removal  of  the  7rf]^(i?123)  contribution,  which  was  shown  to  be 
the  least  significant,  changes  the  spectra  from  a  flat  one,  [7r'^"^(i?1234)]  to  one  that 
has  high  energy  in  the  low  wavenumbers. 
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M,k,5 


Figure  4.12;  One-dimensional  spectra  of  the  total  pressure  compared  with  various 
combinations  of  source  terms  and  regions.  All  regions  and  source  terms; 

.  C4 - C17;  A  C27. 


Figure  4.13;  One-dimensional  spanwise  spectra  of  the  TT  pressure  for  regions  1 
through  3  with  various  source  terms  removed.  All  source  terms;  A  (77123) 
removed;  o  7r^^(i?123)  removed. 


76 


4.5  Summary 


The  MS  and  TT  terms  are  of  the  same  order  of  magnitude  throughout  the  wavenum¬ 
ber  range;  both  have  the  same  shape  as  the  total  pressure.  The  medium  and  high 
wavenumber  {kxS,  kzS  >  5)  MS  partial  pressure  is  generated  in  the  viscous  shear- 
layer  and  buffer  layer,  whereas  the  lower  wavenumber  {kxS,  kz5  <  5)  MS  partial 
pressure  is  generated  in  the  buffer  region  and  logarithmic  region.  The  contribution 
of  the  viscous  shear-layer  to  the  low  wavenumber  MS  partial  pressure  is  insignif¬ 
icant  while  the  upper  channel  contribution  to  the  low  wavenumber  MS  partial 
pressure  is  small. 

The  medium  and  high  wavenumber  TT  partial  pressTires  are  generated  in  the 
viscous  shear-layer  and  in  the  buffer  layer.  In  the  low  wavenumbers  the  TT  par¬ 
tial  pressure  is  due  to  the  buffer  region  and  logarithmic  region  with  a  significant 
contribution  from  the  viscous  shear-layer.  In  the  lowest  wavenumbers  {kx5  <  1), 
all  regions  of  the  channel  are  important.  The  dominant  TT  partial  pressures  in 
the  medium  waveniimber  range  are  due  to  and  from  the  buffer  region. 
In  the  highest  streamwise  wavenumbers  is  about  the  same  magnitude  as 
and  TXi^. 

The  very  highest  spanwise  wavenumber  {k^S  >  60)  partial  pressures  are  domi¬ 
nated  by  TT^^,  7X22  and  txJ^  from  the  viscous  shear-layer. 

The  individual  TT  terms  and  the  TT  partial  pressures  for  various  regions  are 
very  high  in  the  low  wavenumbers.  This  is  in  contrast  to  the  total  TT  spectra 
which  is  flat  in  the  low  wavenumbers.  It  appears  that  the  the  low  wavenumber 
sources  are  coherent  across  the  channel. 
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Chapter  5 


Characteristics  of  Source  Terms 


In  the  previous  Chapter  the  characteristics  of  the  partial  pressures,  due  to  the 
individual  terms  from  four  regions  of  the  channel,  were  examined.  From  that  in¬ 
vestigation,  the  dominant  terms  and  regions  which  contribute  to  the  wall  pressure 
were  discovered.  In  this  chapter,  the  source  terms  themselves  are  studied,  par¬ 
ticularly  those  related  to  the  dominant  partial  pressures.  The  examination  will 
take  place  primarily  in  one-dimensional  spectral  space  at  certain  characteristic 
wall-normal  locations,  which  will  give  indications  as  to  the  most  important  length 
scales  of  the  source  terms.  These  length  scales  will  then  be  related  to  current 
notions  of  boundary-layer  turbulence  structures. 

The  source  terms,  Tij{x,y,  z),  were  computed  in  three-dimensional  physical 
space  using  the  2/3  rule  to  eliminate  aliasing  errors.  Nine  terms  were  com¬ 
puted;  the  six  TT  terms,  their  sum,  the  MS  term  and  the  total.  They  were 
then  Fourier-transformed  to  yield  their  complex  coefficients  in  planar  wavenumber 
space,  The  two-dimensional  spectrum  of  the  source  terms  is  given  by 

^ij{K,y)  =f*{K,y)fij{K,y)  (no  summation)  (5.0.1) 
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(b) 

Figure  5.1:  Source  term  spectra  at  =  5  (a)  streamwise  (b)  spanwise;  A  : 

.  nrpTT.  .  r\rpTT .  w  .  'T'TT.  _  .  r)rpTT .  .  r-pTT .  .  r-pMS. 

.  •  ^-^12  )  .  •  ^-tl3  )  ^  •  -‘22  )  •  ^-‘23  )  •  -‘SS  5  •  i 


_ .  ij^TT .  j-iTot 

The  one-dimensional  spectra  (f)ij{kx,y)  and  4>ij{kz,y)  were  obtained  by  integra¬ 
tions,  similar  to  (4.1.7)  and  (4.1.8)  and  normalized  by  The  one- dimensional 

spectra  of  all  the  source  terms,  at  representative  wall-normal  locations,  are  shown 
in  Figures  5.1  through  5.5.  The  locations  are  in  the  viscous  shear-layer  (y+  =  5); 
in  the  buffer  layer  {y^  =  12  and  21),  which  correspond  to  the  peaks  of  the  m.s. 
distributions  of  the  MS  and  total  source  terms,  respectively;  and  two  locations  in 
the  logarithmic  region,  one  in  the  inner  part,  y^  =  49,  and  one  near  the  channel 
centerline,  y^  =  145.  We  will  now  discuss  the  physics  of  the  MS  term,  followed  by 
a  discussion  of  the  TT  terms. 


5.1  MS  sources 
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Figure  5.2:  Source  term  spectra  at 
legend  as  Figure  5.1. 


(b) 

12  (a)  streamwise  (b)  spanwise.  Same 


Figure  5.3;  Source  term  spectra  at  y'^  =  21  (a)  streamwise  (b)  spanwise.  A  ; 
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Figure  5.4:  Source  term  spectra  at 
legend  as  Figure  5.3. 


(b) 

49  (a)  streamwise  (b)  spanwise.  Same 


Figure  5.5;  Source  term  spectra  at  =  145  (a)  streamwise  (b)  spanwise.  A  ; 


Figure  5.6;  Wall-normal  distribution  of  the  components  of  the  MS  source  term. 
-|-  root-mean-square  values;  o  dv/dx]^—  [dU / dy){dv / dx).  All  values 

have  been  normalized  by  Ur  and  5. 

In  the  previous  Chapter,  the  MS  term  was  shown  to  be  greater  than  the  TT  term 
over  the  entire  streamwise  and  most  of  the  spanwise  wavenumber  ranges.  It  was 
also  shown  that  the  largest  contribution  to  the  MS  partial  pressure  came  from  the 
buffer  layer.  The  MS  term,  shown  in  Figure  5.6,  is  twice  the  product  of  dU/dy  and 
dv /dx.  While  dU/dy  is  maximum  at  the  wall,  dv/dx  is  zero  at  the  wall,  and  peaks 
at  ~  60.  Since  their  product  peaks  at  =  12,  the  most  important  region  for 
MS  partial  pressure  generation  is  expected  to  be  the  buffer  layer. 


Figures  5.1  through  5.5  show  that  the  MS  and  TT  terms  have  fundamentally 
different  behavior  as  kx  0;  the  TT  terms  approach  a  constant  while  the  MS 
term  is  proportional  to  k^.  The  mathematical  reason  for  this  can  be  seen  through 
the  spectral  representation  of  the  MS  term; 


(f)Ms{kx,y)  =  4 


dv  dv 
dx  dx 


=«(f 


(5,1,2) 


which  tends  to  0  as  ^  0  because  both  dU/dy  and  |F*F|  are  well-behaved  func¬ 
tions  as  /ct  — ^  0. 
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(a)  (b) 

Figure  5.7;  One-dimensional  spectra  of  the  MS  source  term.  -|-  y"'"  =  5;  - 

=  12;  A  y~^  —  21;  o  y+  =  49;  a  y+  =  145.  (a)  Streamwise;  (b)  spanwise. 

Comparing  the  spectra  from  the  various  wall-normal  locations  in  Figure  5.7 
shows  that,  consistent  with  Figure  5.6,  the  highest  levels  of  the  MS  source  spectra 
are  found  at  =  12  and  =  21.  This  is  consistent  with  the  finding  that  the 
MS  partial  pressure  from  the  buffer  layer  is  dominant.  The  spectra  at  y”'"  =  5 
and  =  12  have  a  broad  peak  centered  at  approximately  kxS  =  8,  kzS  =  13.5, 
indicative  of  structures  that  are  about  140  and  60  wall  units  in  the  streamwise  and 
spanwise  directions,  respectively. 

Several  observations  suggest  that  near-wall  shear  layers  may  be  the  principal 
source  of  the  MS  pressure;  among  them  are  the  fact  that  near-wall  shear  layers  and 
spanwise  vortices  have  large  gradients  in  the  x  —  y  plane,  which  means  that  dv/dx 
should  be  large  there.  It  has  been  shown  that  the  near-wall  region  is  well-populated 
with  near-wall  shear  layers  but  not  by  spanwise  vortices  [Robinson  et  al.  (1989)]. 
Also,  experimental  and  computational  studies  have  shown  that  they  generate  high 
amplitude  WPFs  [Thomas  and  Bull  (1983),  Johansson  et  al  (1987a,b),  Johans¬ 
son,  Alfredsson,  and  Kim  (1988),  Wilczynski  (1992)].  Finally,  quasi-streamwise 
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vortices,  the  dominant  vortex  structure  in  the  near-wall  region,  do  not  directly 
generate  MS  partial  pressure. 

Near-wall  shear  layers  are  turbulent  structures  formed  at  the  interface  between 
a  downstream  region  of  slower-moving  fluid  and  an  upstream  region  of  faster- 
moving  fluid.  They  are  thought  to  form  when  faster-moving  outer  fluid  impinges 
on  slower  moving  fluid  that  is  being  ejected  by  either  the  legs  of  a  quasi-streamwise 
vortex  or  the  spanwise  head  of  the  horseshoe  vortex.  Johansson  et  al  (1987a) 
showed  that  they  exist  primarily  in  the  near  wall  region  {y^  <  100);  that  in  the 
viscous  shear-layer,  near-wall  shear  layers  form  a  small  (~  7°)  angle  to  the  wall, 
whereas  further  from  the  wall  assume  a  shallow  (~  20*^)  angle  [Figure  5.8(a)]. 
Figure  5.8  shows  that  the  near-wall  shear  layer  has  an  “ejection”  region  (Q2)* 
downstream  of  the  shear  layer  and  a  “sweep”  region  (Q4)  upstream.  Figure  5.8(b) 
shows  that  there  are  three  regions  of  dvfdx:  a  region  in  the  center  of  the  shear 
layer,  where  there  is  a  sharp  transition  from  a  region  of  positive  v,  downstream 
and  negative  v,  upstream;  regions  upstream  and  downstream  of  the  shear  layer 
(?/+  ±125).  A  positive  pressure  event,  shown  in  the  x  —  z  plane.  Figure  5.9, 

is  located  beneath  the  center  of  the  shear  layer.  There  is  very  close  correlation 
between  the  positive  dv/dx  and  positive  pressure,  as  shown  in  Figure  5.10.  The 
relationship  between  and  partial  pressure  at  the  wall  is 

pMS(^y  =  -l,V,K)  -  ^|^;^exp[-A:(l  ±r/)]  (5.1.3) 

where  all  factors  except  dv/dx  have  to  be  positive,  meaning  that  the  sign  of  the 
MS  partial  pressure  is  the  same  as  dv/dx.  Figure  5.10  shows  that  there  is  a  close 
correspondence  between  regions  of  positive  dv/dx  and  positive  pressure.  Similar 

*The  Qi  denote  quadrants  of  the  u  —  v  plane:  Q2  denotes  negative  u,  positive  v;  Q4  denotes 
positive  u,  negative  v.  See  Wallace,  Eckehnann  and  Brodkey  (1972). 
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(a) 


-SOO  -200  -100  0  too  200  300 

(b) 

Figure  5.8:  Shear  layer  structure  as  represented  by  velocity  contours  in  the  x  —  y 
plane.  Horizontal  and  vertical  coordinates  are  in  wall  units.  Velocities  normalized 

by  their  local  RMS  values;  -  :  positive;  .  :  negative,  (a)  u,  contour 

increment  0.5  (b)  u,  contour  increment  0.1.  From  Johansson,  Alfredsson,  and 
Kim  (1988);  published  with  permission. 

results,  in  the  time  domain,  have  been  shown  by  Wilczynski  (1992). 

The  question  remains,  do  the  dominant  length  scales  indicated  by  the 
source  spectra  compare  with  those  for  near-wall  shear  layers?  As  seen  in  Fig¬ 
ure  5.10,  the  streamwise  length  scale  of  the  central  (—60  <  x"*"  <  80)  u-pattern 
is  approximately  140  wall  units.  The  MS  source  spectra  also  has  a  peak  corre¬ 
sponding  to  a  wavelength  of  about  140  wall  units.  Also,  the  MS  partial  pressure 
from  the  buffer  layer,  shown  in  Figure  4.5(a),  has  a  peak  at  kx5  =  3.333  which 
corresponds  to  a  streamwise  length  scale  of  about  380  viscous  units,  close  to  the 
streamwise  length  of  the  entire  shear-layer  structure  and  pressure  pattern  shown 
in  Figure  5.10. 

It  has  been  shown  by  Wilczynski  (1992)  that  the  upstream  (x+  <  —60)  and 
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Figure  5.9;  Wall-pressure  due  to  shear  layer  in  the  x  —  z  plane.  Pressure  normalized 

by  Tyj\  contours  every  0.1;  -  :  positive;  .  :  negative.  From  Johansson, 

Alfredsson,  and  Kim  (1988);  published  with  permission. 


Figure  5.10;  Relationship  between  vertical  velocity  and  wall  pressure  in  centerplane 
of  near- wall  shear  layer  (|/+  =  15).  Horizontal  coordinate  are  in  viscous  units  with 

=  0  detection  location.  ;  v  normalized  by  local  r.m.s.  values;  .  ; 

wall  pressure  normalized  by  r^-  Data  extracted  from  Johansson,  Alfredsson,  and 
Kim  (1988). 
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downstream  {x'^  >  80)  regions  may  have  much  smaller  dvjdx  than  in  the  center 
of  the  shear  layer.  Furthermore,  the  regions,  upstream  and  downstream  of  the 
center  may  contribute  both  to  the  higher  wavenumber  sources,  appropriate  to 
their  local  length  scales  (i.e.,  approximately  100  wall  units),  as  well  as  to  lower 
wavenumber  sources,  which  reflect  the  entire  400  wall  unit  length.  When  the 
source  is  multiplied  by  the  Green’s  function  to  obtain  the  wall  pressure,  the  low 
wavenumber  contributions,  due  to  the  entire  structure  are  magnified  and  show  up 
as  a  lower  wavenumber  peak  in  the  partial  pressure  spectra,  whereas  the  higher 
wavenumbers  are  attenuated.  This  may  explain  why  the  source  spectra  has  a  peak 
that  reflects  only  the  central  part  of  the  dv/dx  signature,  while  the  partial  pressure 
spectra  reflects  the  entire  length  of  the  signature. 

The  spanwise  length  scale  for  the  near-wall  shear  layers  is  about  70  —  80  wall 
units  [Johansson,  Alfredsson,  and  Kim  (1988)].  This  corresponds  well  with  the 
peak  of  the  MS  source  spectra  (~  60  wall  units). 

Overall,  the  length  scales  of  the  streamwise  and  spanwise  wavenumber  spectra 
correspond  very  well  with  the  length  scales  from  the  conditionally  averaged  near¬ 
wall  shear  layers,  providing  more  evidence  that  the  MS  pressure  is  caused  by 
near- wall  shear  layers. 


5.2  TT  sources 

In  this  section  we  attempt  to  make  connections  between  the  dominant  TT  source 
terms  and  coherent  turbulent  structures. 
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Figure  5.11:  Conceptual  drawing  of  a  horseshoe  vortex. 

Quasi-streamwise  vortices 

We  make  the  conjecture  that  the  dominant  TT  terms  are  generated  by  vortical 
structures  such  as  horseshoe  vortices.  As  shown  in  Figure  5.11,  such  a  model 
involves  a  spanwise  vortex  “head”  with  a  tilted  “neck”  and  trailing  “legs”  that 
descend  into  the  buffer  region.  The  vortices  that  make  up  the  legs,  which  have 
a  streamwise  orientation,  and  the  necks,  which  have  a  wall-normal  component 
as  well,  are  called  “quasi-streamwise”.  It  has  been  shown  by  Robinson  (1991b) 
that  the  centers  of  quasi-streamwise  vortices  occur  predominantly  in  the  region 
20  <  <  70,  ~  30  being  their  average  location;  their  diameters  are  in 

the  range  15  <  <  50,  and  their  streamwise  scale  is  “an  order  of  magnitude 

shorter  than  the  longest  sublayer  streaks”  (which  are  approximately  1000  wall 
units)  [Robinson  (1991b)]. 

It  was  found  that  from  the  buffer  layer  is  the  largest  term  over  almost 
the  entire  kx  and  kz  ranges.  The  corresponding  source,  ^ 
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can  be  related  to 


quasi-streamwise  vortices,  which  have  regions  where  the  product  of  dw/dy  and 
dv /dz  is  large;  for  a  circular  vortex,  this  will  occur  in  four  regions  centered  around 
the  quadrant  bisectors,  as  shown  in  Figure  5.12a.  For  instance,  if  a  vortex  is 
centered  at  ~  30  and  has  a  diameter  25,  the  regions  of  maximum  T2^' 

will  occur  at  y^  ~  21  and  y'^  ~  38.  The  first  location  is,  in  fact,  fairly  close  to 
the  peak  of  the  mean-square  distribution  shown  in  Figure  5.13.  The  spectra 
of  at  y^  =  21,  shown  in  Figure  5.3,  shows  that  (kx)  is  dominated  by  low 
wavenumbers  while  '1^2^{kz)  has  a  peak  at  about  kzS  =  25  corresponding  to  long 
streamwise  structures  with  a  spanwise  scale  of  45  wall  units.  This  is  a  reasonable 
length  scale  for  generated  by  quasi-streamwise  vortices. 

It  was  shown  that  (i?2)  was  slightly  larger  than  (i?2)  in  the  interme¬ 
diate  wavenumber  range.  The  term  would  be  large  in  a  vortex  oriented  at 
some  angle  to  the  wall,  such  as  occurs  as  the  quasi-streamwise  vortices  lifts  away 
from  the  wall,  forming  the  “neck”  region  of  the  horseshoe  vortex,  as  shown  in 
Figure  5.12(b).  As  the  angle  of  inclination  increases,  replaces  T'2^^  as  the  pre¬ 
dominant  contribution  of  the  vortex  to  the  soirrce  terms.  Thus,  the  region  of  T'f^ 
generation  should  be  further  from  the  wall  than  for  ,  a  conjecture  confirmed 
by  Figure  5.13,  which  shows  that  the  peak  of  the  mean-square  distribution  of 
is  at  about  y^  =  25.  The  streamwise  wavenumber  spectra  of  at  =  12, 
Figure  5.3,  has  a  peak  at  kxS  ~  10  which  corresponds  to  a  streamwise  wavelength 
of  about  110  wall  units.  This  is  a  reasonable  wavelength  for  a  pressure  footprint 
generated  by  the  neck  region  of  the  horseshoe  vortex  (  e.g.,  a  neck  region  extending 
from  y^  =  20  out  to  y^  =  120,  at  a  45°  angle,  would  generate  a  TTia  footprint 
approximately  100  wall  units  long). 
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quasi-streamwise  vortex 


logarithmic  region 
buffer  region 


viscous  shear  layer 


Figure  5.12;  Schematic  showing  the  relationship  between  horseshoe  vortices  and 
source  terms,  (a)  End  view  showing  showing  regions  of  high  T'£^ ,  and 

(b)  side  view  showing  regions  of  high  and  ■ 

Near- wall  shear  layers  and  spanwise  vortices 


In  Chapter  4,  it  was  shown  that  from  the  buffer  layer  is  one  of  the  dominant 
terms  for  20  <  <  30.  is  the  product  of  du/dy  and  dv/dx  which  may  be 

large  in  near-wall  shear  layers  or  spanwise  vortices.  The  spectra  of  at  various 
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Figure  5.13;  Profiles  of  r.m.s.  TT  and  MS  source  terms.  A  :  .  : 

.  r^rpTT .  V,  .  ^^TT.  .  r^rpTT .  .  t^TT.  _  .  mMS 

.  •  ^-^13  )  ^  ■  -‘-22  I  ■  ^-‘-23  1  •  -‘33  1 

wall-normal  locations,  shown  in  Figure  5.14,  highlight  two  different  behaviors  of 
0”: 


1.  At  the  inner  locations,  =  5, 12  and  21,  the  streamwise  spectra  has  a  peak 
at  about  =  10,  corresponding  to  a  AJ  ~  110.  The  spanwise  spectra 
at  =  5  has  a  peak  at  kzS  =  20,  corresponding  to  a  ~  60;  as  y^  is 
increased  the  peak  in  the  spanwise  spectra  disappears,  and  the  energy  slowly 
moves  to  the  lower  wavenumbers.  The  length  scales  for  can  be  inferred 
from  the  contour  plots  of  uv  from  Johansson,  Alfredsson,  and  Kim  (1988). 
They  show  that  near-wall  shear  layers  have  a  region  of  high  ut)  beneath  the 
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(a)  (b) 

Figure  5.14:  One-dimensional  spectra  of  the  T[.^  source  term.  -|-  =  5;  - 

=  12;  A  ?/+  =  21;  o  y+  =  49;  □  y+  =  145.  (a)  Streamwise;  (b)  spanwise.  Note 
exaggerated  vertical  axis  (compared  with  Figures  5.1  through  5.5  and  5.7). 

ejection  side  of  the  shear  layer  from  which  it  can  be  estimated  that  the  region 

of  high  T'f^  has  streamwise  and  spanwise  length  scales  on  the  order  of  140 

and  60  wall  units,  respectively.  These  wavelengths  are  consistent  with  the 

notion  that  near-wall  shear  layers  generate  However,  at  e.^.,  =  12, 

is  smaller  than  by  about  8  dB;  the  reason  for  this  may  be  that 

du/dy  provides  much  weaker  amplihcation  for  dv/dx  than  does  dU/dy. 

2.  In  Figure  5.14  the  character  of  the  spectra  changes  for  the  outer  two  locations 
(y+  =  49  and  145).  It  can  be  seen  that  the  streamwise  spectra  has  a  slightly 
more  accentuated  peak  which  is  centered  about  kxS  =  12.  This  represents  a 
streamwise  wavelength  of  about  75  wall  units.  The  spanwise  spectra  shifts 
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to  lower  wavenumbers  as  distance  from  the  wall  increases;  this  means  that 
the  dominant  spanwise  wavelengths  are  very  long.  This  behavior  may  be  due 
to  spanwise  vortices.  The  statistics  from  Robinson  et  al.  (1989)  show  that 
most  of  the  spanwise  vortices  occur  for  80  <?/■*■<  180,  although  there  are 
a  significant  number  closer  to  wall  {e.g.,  at  =  50).  The  average  spanwise 
vortex  diameter  at  y'^  =  150  is  about  60  wall  units,  which  compares  well 
with  the  streamwise  wavelength  of  at  y~^  =  49  and  145.  However,  since 
these  sources  have  relatively  low  energy  levels  (compared  with  the  MS  term) 
and  are  far  from  the  wall;  they  are  highly  attenuated  by  the  Green’s  function 
and  generate  relatively  small  partial  pressures. 

Impinging  sweeps 

As  shown  in  Figure  4.11,  the  smallest  scales  of  pressure  are  dominated  by 
and  7^22  due  to  sources  in  the  viscous  shear-layer.  and  are  comprised  of 
the  wall-normal  derivative  of  the  wall-normal  velocity  and  the  spanwise  derivative 
of  the  spanwise  velocity,  respectively.  Such  terms  may  be  related  to  the  impinge¬ 
ment  of  relatively  high-velocity  fluid  on  the  wall  that  occurs  on  the  downward 
side  of  a  quasi-streamwise  vortices  as  shown  in  Figure  5.12(a).  Impinging  flows 
would  create  a  high  wall-normal  gradient  of  the  wall-normal  velocity,  due  to  the 
deceleration  of  the  flow  as  it  impinges  on  the  wall,  and  a  high  spanwise  gradient 
as  the  flow  spreads  out,  parallel  to  the  wall.  could  also  be  generated  as  the 
flow  makes  the  90°  bend. 

The  spectra  of  and  at  y^  =  5,  shown  in  Figure  5.15  show 

that  the  source  terms  are  dominated  by  low  k^,  indicative  of  long  streamwise 
structures.  In  the  spanwise  spectra  both  and  T22'  have  peaks  at  th  50, 
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which  indicates  that  they  have  significant  energy  in  wavelengths  of  about  20  wall 
units;  this  is  a  reasonable  length  scale  for  the  impingement  described  above. 
has  very  little  low  wavenumber  energy,  but  has  a  pronounced  peak  at  th  35, 
which  corresponds  to  a  spanwise  wavelength  of  about  30  wall  units.  The  T2^' 
spectra  at  y"''  =  5  is  about  17  dB  lower  than  the  spectra  for  y"'"  =  21,  shown  in 
Figure  5.3  and  its  peak  has  shifted  from  ^  25  to  k^S  ~  32  which  may  indicate 
that  the  quasi-streamwise  vortices  which  extend  down  into  the  viscous  shear-layer 
have  a  smaller  spanwise  extent  and  may  include  a  contribution  from  secondary 
vortices  formed  during  impingement. 

Robinson  (1991b)  calls  near-wall  regions,  characterized  by  abruptly  diverging 
streamlines,  pockets;  their  shape  is  roughly  circular  in  the  x  —  z  plane  with  span- 
wise  dimensions  of  50  to  100  wall  units;  they  can  be  regions  where  small-scale 
secondary  vortices  form.  It  was  not  clear  to  Robinson  (1991b)  what  the  genera¬ 
tion  mechanism  for  pockets  was,  although  sweeps  were  the  most  likely  candidate 
for  their  generation.  Our  data  indicate  that,  in  particular,  and  are  much 
longer  in  the  streamwise  direction  and  narrower  in  the  spanwise  direction  than 
indicated  by  Robinson’s  (1991b)  description  of  pockets  and  thus,  they  may  not  be 
generated  by  the  same  mechanisms. 


5.3  Summary 

In  summary,  the  MS  source  spectra  strongly  indicates  that  the  MS  pressure  is 
generated  by  near-wall  shear  layers.  Both  the  dominant  length  scales  from  the 
the  MS  spectra  and  the  MS  partial  pressure  compare  very  well  with  data  from 
conditionally  sampled  near-wall  shear  layers.  Since  previous  studies  have  linked 
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Figure  5.15:  One-dimensional  spectra  of  source  terms  from  the  viscous  shear-layer, 

=  5.  -  ^  '^W'l  °  ■  (^)  Streamwise  (b)  spanwise. 

near-wall  shear  layers  to  the  occurence  of  high-amplitude  wall  pressure  events,  it 
can  be  hypothesized  that  the  MS  partial  pressure  is  primarily  generated  by  these 
events. 

The  dominant  TT  terms  occur  in  the  buffer  layer.  The  terms  which  contribute 
most  to  the  wall  pressure  are  and  which  occur  in  quasi-streamwise 

vortices  that  are  parallel  to,  or  at  an  angle  to  the  wall.  T'r2^,  which  was  shown  in 
the  previous  chapter  to  make  only  minor  contributions  to  the  pressure,  is  due  to 
near-wall  shear  layers  in  the  buffer  layer,  and  spanwise  vortices,  in  the  logarithmic 
region.  In  both  cases,  however,  it  is  not  a  very  efficient  generator  of  wall  pressure. 

In  the  viscous  shear-layer  and  are  important  and  appear  to  be 

generated  by  the  downward  side  of  quasi-streamwise  vortices  which  impinge  on  the 
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wall.  These  account  for  the  highest  spanwise  wavenumber  partial  pressures. 


96 


Chapter  6 


Some  aspects  of  MS  pressure  modelling 


One  practical  application  of  the  knowledge  gained  from  the  study  of  WPFs  is  in  the 
field  of  flow-induced  radiated  noise;  noise  prediction  is  of  importance  in  the  design 
and  analysis  of  turbomachinery,  where  the  blade  trailing  edges  are  a  primary  noise 
source.  Radiated  noise  is  related  to  the  wall  pressure  through  an  Helmholtz  equa¬ 
tion  with  a  solution  that  is  the  surface  integral  of  the  wall  pressure  spectral  density 
and  the  wall-normal  derivative  of  a  Green’s  function  for  the  surface  [Blake,  Lee 
and  Zawadski  (1998)].  At  the  design  stage,  which  relies  extensively  on  the  solution 
of  the  Reynolds-averaged  Navier-Stokes  (RANS)  equations  to  obtain  the  velocity 
fields,  to  obtain  the  radiated  noise,  some  model  must  be  employed  that  uses  only 
the  averaged  quantities  (velocity,  pressure  and  Reynolds  stresses)  that  are  avail¬ 
able  from  RANS  calculations.  Such  models  were  developed  by  Blake  (1971,1986), 
Chase  (1980)  and  Zawadski  et  al.  (1996),  using  experimental  data,  theoretical  con¬ 
siderations  and  scaling  arguments.  The  database  generated  in  this  study  allows  us 
to  compute  term-by-term  the  assumptions  required  by  these  models,  and  evaluate 
them  directly. 

An  expression  for  the  MS  pressure  spectrum  for  a  turbulent  one-dimensional 
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mean  flow  was  developed  first  by  Lilley  (1960)  and  modified  by  Blake  (1971,1986), 
Zawadski  et  al.  (1996)  and  Blake  et  al.  (1998).  In  this  formulation  the  spectrum  is 
given  in  terms  of  the  double  wall-normal  integral  of  the  product  of  the  mean  shear 
{dU I dy)  and  the  wall-normal  turbulence  intensities,  with  additional  factors  such  as 
the  one-dimensional  streamwise  and  spanwise  wavenumber  spectra  and  wall-normal 
correlation  of  the  wall-normal  velocity.  The  mean-shear  can  be  obtained  directly 
from  RANS  computations  with  the  wall-normal  turbulence  intensities  estimated 
from  the  RANS  turbulence  model.  The  shortcomings  of  this  model  are:  first, 
it  assumes  one- dimensional  mean  flow;  second,  only  the  MS  term  is  modelled, 
whereas  the  TT  term  is  neglected.  Regarding  the  first  shortcoming,  the  MS  term 
is  the  expansion  of  for  a  three-dimensional  mean  flow  this  results  in  nine 

terms.  As  will  be  seen,  there  is  no  accurate  model  even  for  the  single  term  that 
appears  for  one-dimensional  mean  flows.  Thus,  part  of  a  modelling  effort  would 
be  to  determine  how  to  apply  a  one-dimensional  model  in  a  three-dimensional 
flow.  Regarding  the  second  shortcoming,  it  was  shown  in  Chapter  4  that  the  one¬ 
dimensional  MS  pressure  spectra  is  slightly  larger  than  the  TT  pressure  at  all  but 
the  highest  spanwise  wavenumbers,  and  its  spectra  has  the  same  shape  as  the  total 
pressure,  and  is  lower  by  a  few  decibels.  Thus,  the  MS  pressure  may  provide  a 
reasonable  approximation  to  the  total  pressure. 

In  this  chapter,  we  review  the  formulations  of  Blake  (1986),  Chase  (1980)  and 
Zawadski  et  al.  (1996),  restate  their  formulations  for  channel  flow,  and  evaluate 
their  various  levels  of  approximation.  Then,  we  shall  study  the  one-dimensional 
streamwise  and  spanwise  spectra  and  broadband  wall-normal  correlation  coeffi¬ 
cient  of  V,  which  are  factors  in  the  formulations  that  can  only  be  modelled  with 
knowledge  of  turbulent  velocity  fields,  and  are  conveniently  computed  from  the 
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DNS  database. 


6.1  MS  model  formulation 


The  solution  to  the  Poisson  equation  for  fluctuating  pressure  for  channel  flow  was 
derived  in  Section  2.6.  Ignoring  the  boundary  terms  (in  Section  4.2  it  was  shown 
that  the  Stokes  pressure  was  very  small  compared  with  the  total  pressure,  except  at 
the  highest  spanwise  wavenumbers),  the  solution  is  given  in  wavenumber-frequency 
space  by 


p(y,  K,uj)  =  -  l'f‘-‘{v,K,u,)g{y,K,ui)dy. 

where  the  Green’s  function  evaluated  at  the  near-wall  is  given  by 


9{K,  y) 


cosh[iT(l  —  y)] 
Ksm\v(2K) 


(6.1.1) 


(6.1.2) 


The  only  differences  between  the  channel  flow  and  boundary  layer  solutions  is  the 
Green’s  function,  which  for  a  boundary- layer  is  and  the  limits  of  integration, 
[0,oo]. 


The  spectral  density  (amplitude)  of  the  wall  pressure  is  given  by  multiplication 
of  the  complex  spectra  (6.1.1)  by  its  complex  conjugate 

\p*{K,uj)  p{K',u')\ 


n'°*(K,a;)  = 


(6.1.3) 


d  (K  -  K')  5{uj-  u')  ■ 

The  6  functions  indicate  that  the  function  need  only  be  evaluated  at  K  =  K'  and 
to  =  uj'.  The  spectral  density  for  the  wall  pressure  in  terms  of  (6.1.1)  is  given  by 

1 

dy  g{y)  \f^^*{y)  +  f^^^*{y)\  J  dy'g[y')  [f^'^[y')  +  f^^^{y') 

(6.1.4) 


-1 
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It  is  understood  that  the  arguments  in  the  right-hand-side  of  (6.1.4)  are  actually 
(y,  K,  uj)  or  {y',  K,a;).  The  wall-normal  integration  in  (6.1.1)  must  be  over  separate 
dummy  integration  variables,  leading  to  a  double  integration  in  the  wall-normal 
direction.  One  could  also  evaluate  p*(0,  K,a;)  and  p(0,K',a;')  separately,  then 
take  their  amplitude,  thus  avoiding  the  complications  involved  with  the  double 
wall-normal  integration.  Such  an  approach,  however,  does  not  lead  to  a  form 
involving  the  product  of  the  wall-normal  correlations  and  turbulence  intensities, 
which  make  the  MS  term  at  all  possible  to  model.  Multiplication  of  the  terms  in 
(6.1.4)  gives 

1  1 

n*°*(K,  u;)=  I  dy  I  dy'  G{y,  y')  [f  ^^*(y)f^^(yO  +  f^^^*{y)f^^^y')  +  f^\ 

-1  -1 

(6.1.5) 

where  G{y,y')  =  g{y)g{y')  and  represents  the  cross  terms  given  by 

f^{y,y')  =  f^^*{y)f^^^{y')  +  f^^{y')f^^*{y).  (6.1.6) 
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which  is  an  exact  representation  of  the  MS  spectral  density  and  still  contains 
the  turbulent  quantity  v.  The  challenge  for  various  modellers  [e.g..  Chase  (1980), 
Zawadski  et  al.  (1996),  Blake  et  al.  (1998)]  has  been  to  model  v*{y^  K,  uj)v{i/,  K,  uj) 
using  empirical  knowledge  of  turbulence. 

Modelling  the  vertical  velocity 

Blake  (1986)  defines  <h22(|/,  ?/^K  ,u)),  the  cross  spectral  density  of  the  normal  com¬ 
ponent  of  fluctuating  velocity  between  wall  distances  y  and  y',  to  be 

<i>22iy,y',i^,uj)  =  \v*{y,K,u;)v{y',K,u))\ .  (6.1.11) 

If  the  magnitude  operation  in  (6.1.10)  is  brought  inside  the  integral,  then  the 
MS  wall  pressure  spectral  density  (6.1.10)  can  be  expressed  in  terms  of  the  cross 
spectral  density, 

1  1 

n"®(K.u.)  =  4*^  Jdy  Jdy’G{y,y’]T(y)T(y’)iMy,  v\K,i0)  (6,1,12) 

-1  -1 

The  mathematical  difference  between  forms  (6.1.10)  and  (6.1.12)  is  that  in  (6.1.10) 
the  integrals  of  the  complex  spectra  are  taken,  then  the  magnitude,  whereas,  in 
(6.1.12)  the  magnitude  is  taken  first,  then  the  integrals.  This  would  yield  the 
same  results  if  the  product  of  v(K,y)*v{K.,y')  were  real,  as  is  the  case  at  y'  =  y. 
However,  for  y'  ^  y,  the  product  is  complex  and  form  (6.1.12)  ignores  cross  terms 
that  would  appear  in  the  integration  of  a  complex  function.  The  error  involved  in 
this  approximation  will  be  quantified  later. 

A  further  simplification  is  to  express  the  cross  spectral  density  function  in 
product  form 

<^22(1/, I/',  K, a;)  R22{y,y')^22{y,'i<^,(^)-  (6.1.13) 
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$22(1/,  K,  a;)  is  the  spectral  density  of  wall  normal  velocity  at  a  single  level,  y,  while 


R22{y,  y') 


{v{y)v{y')) 

{v‘^{y')) 


(6.1.14) 


is  the  broadband  wall-normal  correlation  function  (so-called  because  it  is  integrated 
over  all  wavenumbers  and  frequencies),  which  contains  the  relationship  between 
the  velocity  at  y  and  y' .  Furthermore,  $22(1/,  K,  a;)  can  be  approximated  as  the 
product  of  three  one- dimensional  functions 


$22(1/,  K, a;)  ^  -  k.  U,)  {v\y))  •  (6-1.15) 

^vvikx)  and  are  the  one-dimensional,  normalized  wavenumber  spectra  of  v, 

defined  below.  —  k^Uc)  is  the  moving-axis  spectrum  which  is  a  functional 

form  for  the  distribution  of  energy  away  from  the  convective  ridge  (in  the  kx  —  uj 
plane)  for  the  energy  in  the  vertical  velocity.  If  Taylor’s  hypothesis  is  invoked, 
then  00”  becomes  a  delta  function  at  each  value  of  lo 


0m (^  kx  kk c)  —  S{uJ  kx  Uc)‘ 


(6.1.16) 


This  considerably  simplifies  the  integration  of  (6.1.12)  over  all  00  to  obtain  the 
planar  wavenumber  spectra,  since 

00 

J  5{lo  —  kxUc)  du)  —  1.  (6.1.17) 

—00 

Invoking  Taylor’s  hypothesis  does  not  affect  the  results  presented  herein  since  we 
are  only  investigating  the  one-dimensional  spectra,  which  are  obtained  by  integra¬ 
tion  of  the  three-dimensional  spectra  over  all  frequencies  and  one  of  the  wavenum¬ 
ber  directions.  Thus,  since  the  convective  ridge  dominates  the  entire  wavenumber- 
frequency  domain,  the  one-dimensional  spectra  reflect  mostly  the  energy  in  the 
convective  ridge.  If  more  detail  about  the  subconvective  wavenumber  pressures 
are  desired  then  a  more  complex  functional  form  for  0™  needs  to  be  determined. 
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The  functions  (j)^„{kx,y)  and  (j)l^{kz,y)  are  the  normalized  spectra, 


7  /,  x_  4>{kx,y) 

(Pw  (  kx  j  ?/)  — 


Jdkx  4^vv  {kxi  y) 


(6.1.18) 


4>lxi^z)y)  is  normalized  in  a  similar  manner.  By  Parseval’s  theorem, 

^  Lx  OO 

(v^iy))  =  j  j  dxv'^{x,y,z)  =  j  dkx  (pvvikx^y).  (6.1.19) 

^^00  0 

Using  (6.1.15)  in  (6.1.13),  and  placing  the  result  in  (6.1.12)  gives  an  approximate 
form  for  the  MS  pressure; 


1  1 

n^^^(K)  ^4kl  J  dy  j  dy'G{y,y')T{y)T{y')  (^v^{y))  R22{y,y')PxvikxWvvikz)- 

-1  -1 

(6.1.20) 

This  will  be  referred  to  as  the  “3-term”  form.  Note  that  RANS  calculations  would 
yield  {v'^{y))  and  r(y),  but  all  the  other  terms  require  further  modelling. 


Comparison  of  approximate  forms 

Figure  6.1  compares  the  exact  (6.1.7),  CSD  (6.1.12)  and  the  3-term  separable 
(6.1.20)  forms.  Results  are  averaged  from  about  80  realizations  of  DNS5.  The 
differences  between  the  various  forms  are  quite  small  and  occur  primarily  at  the 
lowest  wavenumbers.  The  difference  between  the  exact  and  CSD  forms  shows 
the  effect  of  moving  the  magnitude  operation  from  outside  the  integral  (6.1.7) 
to  inside  (6.1.12).  The  difference  between  the  CSD  and  3-term  forms  show  the 
difference  between  evaluating  the  four-dimensional  <F(y,  y',K)  and  the  product  of 
one-dimensional  spectra. 
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(a)  (b)  M 

Figure  6.1;  Comparison  of  computed  with  three  levels  of  approximation. 

(a)  Streamwise,  (b)  spanwise  spectra.  Normalized  by  r'^S.  -  :  Exact  (6.1.7); 

.  :  CSD  (6.1.12); .  :  3-term  (6.1.20). 

6.2  Characteristics  of  vertical  velocity 

In  this  section  we  study  the  characteristics  of  the  one-dimensional  spectra  and 
broadband  correlations  of  v.  Such  information  is  necessary  for  forms  involving 
evaluation  of  the  spectra  at  characteristic  wall- normal  locations  e.g.,  it  may  be 
necessary  to  model  (see  Zawadski  et  al.  (1996))  the  wall-normal  velocity  spectra 
based  on  a  single,  characteristic  location  y  =  i/c,  *-e., 

n^^‘^(K)  ^  j^^dy  {K,y,y')  R{y,y')4>22{yc,K)4>22{yc,k^)  (t)‘^{y))  . 

(6.2.1) 

Traces  of  022(|/,  kx)  and  (t)22{y,  kz)  at  various  y  locations  are  shown  in  Figure  6.2.  It 
is  immediately  apparent  from  the  k^  spectra  that  022 (y,  k^)  is  practically  constant 
with  y.  The  cut  at  y^  =  145  shows  a  small  divergence  at  the  lowest  wavenumbers, 
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Figure  6.2:  Constant  cuts  of  022(y, K).  (a)  Streamwise  (b)  spanwise.  -  : 

=  5; .  :  2/+  =  12; .  :  =  21; - ;  t/+  =  49;  —  ;  =  145. 

but  for  the  purposes  of  modelling,  this  difference  may  be  insignificant.  On  the 

other  hand  the  spectra  shows  marked  differences  with  y,  although  all  the  curves 

appear  to  have  similar  shapes.  The  peak  of  the  spectra  moves  from  higher  to  lower 

values  of  kz  as  y  is  increased. 

Contour  plots  of  the  spectra  in  and  {kz^y'^)  space,  Figure  6.3,  show 

the  distribution  of  energy  as  a  function  of  wavenumber  and  wall-normal  distance. 
Superimposed  on  top  of  the  plot  are  the  maxima  found  using  a  peak-detection 
algorithm.  As  shown  in  Figure  6.2,  the  peaks  are  very  broad,  so  the  trace  of  the 
maxima  is  only  an  approximation.  In  any  case  it  can  be  seen  that  the  spectra 
are  fairly  constant  in  y,  and  it  might  be  safe  to  assume  that  a  characteristic  of 
022(y,  kx)  can  be  taken  from  almost  anywhere  in  the  boundary  layer.  In  the  viscous 
shear-layer,  kz  spectra  have  a  peak  at  about  kz5  =  12;  the  peak  then  migrates  to 
lower  kz  as  y  increases  away  from  the  viscous  shear-layer. 
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Figure  6.3:  Contour  plot  of  spanwise  wavenumber  spectra,  022(yj^z)-  The  thick 
line  denotes  the  trace  of  the  maxima. 


Figure  6.4:  Peak  kz  values  of  022(?/,  k^),  as  a  function  of  y^.  Labels  show  the  values 
of  n  for  the  y"  behavior  of  the  various  zones. 


Region 

y  min 

y  max 

n 

approximation 

1 

0 

6 

-0.03 

0 

2 

6 

50 

-0.35 

-1/3 

3a 

50 

no 

-1.19 

-1 

3b 

no 

180 

-0.87 

-1 

Table  I;  Wall-normal  regions  for  (f)22{y,  k^)  peak  locations,  n  is  the  exponent  in  the 
relationship  k^S  =  y”. 


Figure  6.5;  4>22{y,kz)  in  the  logarithmic  region  plotted  vs.  kzS.  -  :  =  53; 

. :  y^  =  75; .  :  y+  =  100; - :  y+ =  128;  —  :  y+  =  151; - : 

y+  =  174. 

Figure  6.4(b)  shows  the  location  of  the  peak  of  the  spanwise  wavenumber  spec¬ 
tra.  The  curve  presents  sharp  changes  in  behavior  at  y^  =  6,  50  and  110.  Within 
each  layer,  it  can  be  described  by  a  power-law  behavior,  y”,  with  a  constant  ex¬ 
ponent.  The  exponents  obtained  by  fitting  the  data  are  shown  in  Table  I;  an 
approximate  power-law  behavior  is  also  proposed.  Note  that  in  the  logarithmic 
region  and  beyond,  Morrison  and  Kronauer  (1969)  proposed  a  y^^  behavior  for 
two-dimensional  spectra  of  u,  which  is  consistent  with  the  present  data. 
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Figure  6.6:  022('V)^z)  in  the  logarithmic  region,  plotted  vs.  k^v,  normalized  by 

(a)  (b)  {v^{7]))7].  -  y+  =  53;  .  7j+  =  75;  .  =  100; 

- y+  =  128;  —  y+  =  151; - y+  =  174. 

Figure  6.5  shows  the  spanwise  wavenumber  spectra  from  five  locations  in  the 

logarithmic  region;  the  peak  moves  to  lower  wavenumbers  and  becomes  higher  and 

less  pronounced,  as  y  increases.  The  high  wavenumber  tail  moves  in  a  similar 

fashion  toward  lower  wavenumbers  as  y  increases.  This  indicates  that  the  spectra 

may  scale  like  kzy. 

Morrison  and  Kronauer  (1969)  claimed  that  their  u  spectra  could  be  expressed 
in  product  form 

n(a;+,  iF+|/+)  =  f{K+y+)A{tu+,  kf)  (6.2.2) 

where  ^  denotes  normalization  by  viscous  variables  (cu^  =  uju/ur,  =  Kv/ur). 
f  is  called  the  “wave  intensity  function”  and  A  is  the  “wave  strength.”  This 
implies  that  for  a  range  of  wall-normal  locations,  their  data  collapsed  to  a  single 
curve  when  plotted  vs.  K^y^.  We  showed  previously  that  the  peak  locations 
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of  the  spanwise  wavenumber  spectra  scale  accordingly.  Figure  6.5  shows  that  at 
wavenumbers  higher  than  the  peak,  the  spectra  decrease  with  distance  from  the 
wall.  This  would  lead  one  to  believe  that  the  abscissa,  multiplied  hy  rj  {rj  =  \  +  y 
is  the  distance  from  the  wall),  would  collapse  the  high  wavenumber  spectra  to  a 
single  curve.  Figure  6.6(a)  shows  that  the  logarithmic  region  spectra,  normalized 
by  {v'^{y)),  plotted  versus  k^rj  has  all  the  peaks  lined  up  at  a  single  location, 
kzV  =  3.  While  the  lowest  wavenumbers  collapse  to  a  single  curve,  the  spectral 
levels  at  wavenumbers  higher  than  the  peak  tend  to  increase  with  y.  It  would  seem 
that  normalizing  022  by  y  may  collapse  the  data.  Figure  6.6(b)  shows  that  022(^z)? 
normalized  by  {v^{y))  y  moves  the  curves  for  various  y^  closer  together,  but  that 
in  the  high  wavenumbers  the  spectra  values  increase  with  increasing  y^ . 


Broadband  correlation  function 


In  the  3-term  formulation  (6.1.20),  the  broadband  correlation  function  R22{y,y'), 
characterizes  the  relationship  between  the  vertical  velocities  at  different  wall- 
normal  locations.  Zawadski  et  al.  (1996)  showed  that  for  a  range  of  vertical  lo¬ 
cations  8  <  y"*"  <  60,  R22{y,y'),  normalized  by  the  m.s.  value  at  the  fixed  point, 
(?/')),  collapsed  to  a  single  curve  for  y  <  y\  when  plotted  against  y/y'  {y'  held 
fixed,  while  y  varied).  The  correlations  were  computed  by  the  formula. 


R22{y,y’) 


{v{y)v{y')) 

{v‘^{y')) 


(6.2.3) 


from  60  realizations  of  DNS5  data.  Figure  6.7  shows  that  for  y'  in  the  viscous 
shear-layer  and  buffer  layer  this  results  in  R22{y,  y')  values  greater  than  1.0.  When 
plotted  against  y/y' ^  Figure  6.7(b),  the  curves  collapse  for  y/y'  <  1.0  consistent 
with  the  findings  of  Zawadski  et  al.  (1996).  It  can  be  seen  that  for  y/y'  >  1.0 
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(b)  y/y' 

Figure  6.7;  Broadband  wall-normal  correlations  of  v  normalized  by  {v'^iy')) 

y^'  =  5; .  y+'  =  13; .  y+'  =  21; - y+'  =  49;  —  y+'  -- 

Plotted  against  (a)  y^  (b)  y/y' . 


the  curves  diverge  from  each  other,  with  the  greatest  amount  of  divergence  for  the 
smallest  values  of  y' . 


Normalizing  by  the  r.m.s.  values  from  each  of  the  two  levels 


R22{y,y') 


{v{y)v{y')) 


(6.2.4) 


gives  correlations  with  the  property  R22{y,y')  <  1-0.  The  correlations  are  shown 
in  Figure  6.8  for  a  range  of  y'  values  throughout  the  channel.  Figure  6.8(a)  shows 
that  plotting  the  correlations  plotted  against  y/y'  does  not  collapse  them  to  a 
single  curve.  More  work  is  necessary  to  determine  how  the  correlations  might  be 


modelled. 


Narrowband  correlation  function 


The  narrowband  correlation  functions  provide  wall-normal  information  at  partic¬ 
ular  wavenumbers.  They  were  normalized  by  the  m.s.  value  at  the  fixed  point, 


y 


^22  (y,  y',K) 


{v{y,K)v{y',K)) 

(u2(y',K)) 


(6.2.5) 


Figure  6.9  are  plots  of  the  narrowband  correlations  for  four  wavenumbers  plot¬ 
ted  vs.  y^.  The  four  wavenumber  pairs  {kxS,kz6),  (3,4.5),  (7,7.5),  (16,15)  and 
(30,45),  were  selected  to  cover  a  wide  range  of  wavenumbers.  As  can  be  seen  from 
Figure  6.9,  the  correlations  get  narrower  as  K  increases,  a  result  that  is  not  unex¬ 
pected  since  the  higher  wavenumbers  are  attenuated  by  the  Green’s  function  to  a 
greater  extent  than  the  lower  ones.  As  with  the  broadband  correlations,  the  peak 
location  can  be  collapsed  by  plotting  against  y/y'.  Figure  6.10,  plots  of  of  narrow- 
band  correlations  plotted  vs.  y/y',  show  that,  as  with  the  broadband  correlations, 
for  y/y'  <  1,  there  is  collapse  of  the  data,  particularly  for  77  =  22. 
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(c)  (d)  y" 

Figure  6.9;  Narrowband  wall-normal  correlations  of  v  normalized  by  (u^(j/',K)) 

plotted  vs.  y^.  -  y^'  =  21;  .  y^'  =  53;  .  y'^'  =  152.  At 

wavenumbers  (a)  K  =  5.4,  (b)  K  =  10.3,  (c)  K  =  22,  and  (d)  K  =  54. 

6.3  Summary 

In  this  chapter,  we  compared  the  MS  pressure  computed  using  the  various  levels  of 
approximation  from  the  formulation  of  Blake  (1971).  Results  showed  that  moving 
the  magnitude  operation  from  outside  to  inside  the  wall-normal  integrals,  so  that 
the  cross  spectral  density  function  could  be  used,  increases  the  MS  pressure  by  a 
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1  ■  1 
(c)  y/y'  (d)  yiy' 

Figure  6.10;  Narrowband  wall-normal  correlations  of  v  normalized  by  (u^(y',K)) 

plotted  vs.  y/y'.  -  y^'  =  21;  .  y^'  =  53;  .  y^'  =  152.  At 

wavenumbers  (a)  K  =  5.4,  (b)  K  =  10.3,  (c)  K  =  22,  and  (d)  K  =  54. 

few  dB;  expressing  the  cross  spectral  density  function  in  terms  of  products  further 

increases  the  MS  pressure  by  a  few  dB.  In  both  cases,  the  character  of  the  MS 

pressure  spectra  was  not  altered.  It  can  be  concluded  that  the  formulation  of 

Blake  (1971)  is  not  signihcantly  different  from  the  actual  MS  pressure  and  may  be 

a  good  starting  point  for  further  modelling  efforts. 

We  studied  the  characteristics  of  the  one-dimensional  spectra  of  v  and  broad- 
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band  correlation  functions,  which  are  the  functions  which  must  be  modelled  in  the 
MS  pressure  formulation.  We  found  that  the  streamwise  spectra,  normalized  by 
{v'^{y)),  varies  only  slightly  throughout  the  channel  and  that  a  characteristic  spec¬ 
tra  can  be  taken  from  almost  anywhere.  On  the  other  hand,  the  spanwise  spectra, 
normalized  by  {v'^{y)),  varies  across  the  channel,  but  in  a  way  that  can  potentially 
be  modelled;  the  shapes  of  the  spectra  appear  to  be  similar  and  the  wavenum¬ 
ber  location  of  the  spectra  peak  varies  by  powers  of  y,  in  discrete  regions  of  the 
channel.  An  unsuccessful  attempt  was  made  at  hnding  a  normalization  that  would 
collapse  the  spanwise  wavenumber  spectra.  The  broadband  correlation  functions 
were  shown  for  locations  throughout  the  channel.  When  normalized  by  the  m.s. 
value  at  the  fixed  point,  (u^(y')),  the  correlations  appear  to  collapse  for  y/y'  <  1. 
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Chapter  7 


Concluding  Remarks 

7.1  Summary  and  conclusions 

The  motivation  for  this  investigation  was  to  answer  some  longstanding  questions 
about  the  sources  of  wall  pressure.  Namely,  where  in  the  boundary  layer  are  the 
dominant  sources,  what  are  the  physical  mechanisms  of  these  sources,  and  what  are 
the  relative  magnitudes  of  the  mean-shear  (MS)  and  turbulence-turbulence  (TT) 
interaction  terms.  Using  the  velocity  fields  from  a  direct  numerical  simulation 
(DNS)  of  a  fully  developed,  turbulent,  channel  flow,  the  source  terms  were  com¬ 
puted  in  planar  wavenumber  space.  The  convolution  of  the  Green’s  function  with 
the  source  terms,  over  four  horizontal  layers  of  the  channel,  gave  partial  pressures 
for  each  of  the  terms  and  regions.  Studying  the  partial  pressures  gave  insights 
into  which  source  terms  were  dominant  and  where  in  the  channel  they  were  act¬ 
ing.  Studying  the  spectra  of  the  dominant  source  terms  themselves  showed  that 
the  dominant  terms  are  closely  related  to  generally  accepted  notions  of  coherent 
events  in  the  near-wall  region. 

The  practical  motivation  for  this  research  is  to  model  the  wall  pressure  using  the 
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data  and  information  from  the  turbulent  velocity  fields.  This  is  an  area  for  which 
previous  modelling  attempts  [Blake  (1971),  Chase  (1980)]  lacked  detailed  informa¬ 
tion.  One-dimensional  spectra  and  broadband  correlations,  important  components 
of  the  MS  models,  were  evaluated  and  attempts  were  made  to  model  them. 

The  major  findings  from  this  investigation  are  as  follows: 

1.  The  MS  and  TT  terms  are  of  the  same  order  of  magnitude  throughout  the 
wavenumber  range;  both  have  the  same  shape  as  the  total  pressure. 

2.  The  medium  and  high  wavenumber  >  5)  MS  partial  pressure  is 

generated  in  the  viscous  shear-layer  and  buffer  layer,  whereas  the  lower 
wavenumber  {kxS,  kzS  <  5)  MS  partial  pressure  is  generated  in  the  buffer 
layer  and  logarithmic  region.  The  contribution  of  the  viscous  shear-layer 
to  the  low  wavenumber  MS  partial  pressure  is  insignificant,  and  the  upper 
channel  contribution  to  the  low  wavenumber  MS  partial  pressure  is  small. 
It  is  conjectured  that  the  MS  partial  pressure  is  primarily  due  to  near-wall 
shear  layers,  and  that  high  amplitude  events  commonly  found  beneath  such 
turbulence  structures  will  contribute  primarily  to  the  MS  partial  pressure. 

3.  The  medium  and  high  wavenumber  TT  partial  pressures  are  generated  in 
the  viscous  shear-layer  and  buffer  layer.  At  low  wavenumbers,  the  TT  par¬ 
tial  pressure  is  due  to  the  buffer  and  logarithmic  regions  with  a  significant 
contribution  from  the  viscous  shear-layer  as  well.  At  the  lowest  wavenumbers 
{kxS  <  1),  all  regions  of  the  channel  are  important.  The  dominant  TT  partial 
pressures  in  the  medium  wavenumber  range  are  due  to  and  from  the 
buffer  layer.  These  terms  are  strongly  affected  by  tilted  quasi-streamwise  vor¬ 
tices.  The  very  highest  spanwise  wavenumber  {k^S  >  60)  partial  pressures 
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are  dominated  by  and  from  the  viscous  shear-layer.  These 

partial  pressures  may  be  generated  by  downward  side  of  quasi-streamwise 
vortices  impinging  on  the  wall. 

4.  71^2  was  found  to  be  relatively  large  at  high  streamwise  wavenumbers  and 
may  be  related  to  near-wall  shear  layers,  although  its  contributions  to  the 
wall-pressure  are  small  compared  to  the  MS  term. 

5.  TiJ.]"  also  appears  to  be  related  to  spanwise  vortices  found  in  the  logarithmic 
region.  However,  the  distance  from  the  wall  and  the  concentration  of  energy 
at  relatively  high  wavenumbers  precludes  them  from  being  major  contribu¬ 
tors  to  the  wall  pressure. 

6.  Evaluation  of  the  MS  modelling  formulation  of  Blake  (1971)  shows  that  the 
model  gives  a  reasonable  prediction  of  the  one-dimensional  MS  pressure  spec¬ 
tra  when  exact  values  of  the  vertical  velocity  spectra  and  broadband  corre¬ 
lations  are  used.  In  order  for  such  models  to  be  used  as  general  purpose 
predictive  tools,  however,  the  one-dimensional  spectral  densities  and  wall- 
normal  correlation  function  of  vertical  velocity  need  to  be  modelled.  It  was 
shown  that  the  streamwise  spectra  of  v  is  constant  across  the  channel.  The 
peak  of  the  spanwise  spectra,  however,  shifts  towards  lower  wavenumbers 
with  distance  from  the  wall.  The  location  of  the  peak  varies  as  y^, 

and  y~^  in  the  viscous  shear-layer,  buffer  layer,  and  logarithmic  re¬ 
gion,  respectively.  Further  efforts  are  necessary  to  scale  the  amplitude  of  the 
spanwise  v  spectra.  The  broadband  correlations  were  shown  to  collapse  to  a 
single  curve  for  yjy'  <  1.0  when  plotted  versus  y/y'. 
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7.2  Recommendations  for  further  work 


A  number  of  areas  closely  related  to  the  current  work  could  profitably  be  exploited 
to  clarify  further  the  relationship  between  the  velocity  fields  and  the  wall  pressure 
and  provide  information  for  modelling  in  more  practical  geometries: 

1.  Using  a  channel  flow,  the  wavenumber-frequency  spectra  of  the  dominant 
partial  pressures  could  be  obtained.  This  would  provide  data  regarding  the 
importance  of  the  MS  and  TT  term  contributions  to  the  kx  —  uj  spectrum, 
particularly  with  regard  to  the  subconvective  wavenumber  range.  It  may  also 
help  to  determine  the  cause  of  the  high  energy  in  the  low-wavenumber  TT 
partial  pressures.  Study  of  the  kz  —  uj  space  should  be  undertaken  in  order 
to  compare  with  later  streamwise  varying  flows. 

2.  Conditional  averaging  of  the  velocity  helds  could  be  used  to  validate  our 

hypotheses  regarding  the  relationship  between  the  and  near-wall  shear 
layers,  and  between  and  quasi-streamwise  vortices.  In  particular,  by 

generating  a  time  series  of  MS  and  TT  wall-pressure  fields,  the  hypothesis 
that  near-wall  shear  layers  and  thus,  high-amplitude  wall-pressure  events 
generate  the  MS  pressure  could  be  validated. 

3.  Simulation  of  channel  flow  at  higher  Reynolds  number  would  provide  greater 
separation  (in  viscous  units)  across  the  channel.  Such  a  simulation  could  help 
to  answer  questions  about  the  low-wavenumber  energy  in  the  individual  TT 
terms  and  may  answer  questions  about  the  small  energy  in  partial  pressures 
due  to  the  logarithmic  region  spanwise  vortices. 
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4.  Simulation  of  a  boundary  layer  flow  and  the  evaluation  of  the  partial  pres¬ 
sures  could  yield  information  about  the  low  wavenumber  partial  pressures. 
Such  a  study  could  prove  that  the  high  energy  in  the  low  wavenumbers  of 
the  channel  flow  are  due  to  low- wavenumber  coherence  of  the  opposite  walls 
of  the  channel  flow. 

5.  Simulation  of  flows  with  a  backward  facing  step  or  adverse  pressure  gradient 
would  yield  information  about  the  partial  pressures  and  sources  in  other  than 
equilibrium  flows.  Information  from  such  flows  are  very  important  for  de¬ 
velopment  of  airfoil  trailing  edge  noise  models.  Particular  tasks  would  be  to 
compute  the  partial  pressures  in  ~  uj  space;  determine  the  dominant  par¬ 
tial  pressures  and  their  relationship  to  coherent  events;  validate  the  relative 
importance  of  the  MS  and  TT  partial  pressures;  determine  the  relationship 
between  and  near-wall  shear  layers. 
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